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ABSTRACT 

The developoent and liaited testing of a cooputer control of a 
liquid cooled garoent (LOG) inlet teaperature is described. An 
adaptive (cooputer) acdel of the LCG is used to predict the heat- 
removal rates for various inlet tessera tures. The predictions can be 
used to choose an inlet temperature that falls into a comfort-zone of 
the type that was derived by Kuznetz in earlier studies at NASA. 

An experimental system that contains a microcomputer has been 
constructed. LCG inlet and outlet temperatures and the heat- 
exchanger outlet temperature form the inputs to the computer. The 
cooputer commands a valve to control the LCG inlet temperature* The 
adaptive-model-prediction method of control has proved successful 
during tests on a number of subjects where the inlet temperature was 
automatically chosen by the computer. The program can be implemented 
in a microprocessor of a size that is practical for a life-support 
back-pack. 
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I INTRODUCTION 

A number of techniques have been evaluated for providing a 
control algorithm for automatic control of the inlet temperature to a 
liquid cooling garment (LCG)* The closed loop control system 

will require monitoring of appropriate system (physiological and/or 
LCG) parameters t processing of these pa^'ameters in some manner to 
generate a control signal, adjustment of the LCG heat removal rate in 
response to the control signal, and subsequent monitoring and 
adjustment to insure that the desired response is obtained. 

The heat I'emoval requirement for a given LCG is a function of 
the metabolic rate, the desired thermal comfort conditions, and the 
heat exhange with the environment. However, since it is relatively 
easy to minimize heat exchange with the environment, the metabolic 
heat is the primary concern for a practical system. 

The following physiological or system parameters could be 
monitored and used, with appropriate control algorithms, to provide 
automatic control of thermal comfort while using an LCG. Each will 
be briefly discussed to provide the rationale leading to the 
selection of the parameters and control model developed during the 
present research. 

A. Body Temperature Monitoring 

The body temperatures most suitable for providing an input 
signal to a controller are the core temperature as indicated by a 
rectal or ear canal probe and/or various sxin temperature. However, 
there are problems associated with these for use in an actual 
operational EVA system in addition to some physiological drawbacks. 

Monitoring of body temperatures requires the attachment of 
temperature sensors to the subject. This Is undesirable for a system 
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which will be used for hundreds of manhours on each mission. In 
addition, the time constant for the core temperature to respond to a 

•ago 

sudden change in metabolic rate is on the order of 10 minutes'^’ ’ so 

that the control algorithm must anticipate the body’s response to 

changes in inlet temperature. Skin temperature(s) is not a simple 

function of metabolic rate*^ and thus would be difficult to use by 

itself as the only input to eui LOG controller. A combination of core 

0 

and surface (skin) temperature could probably be made to work well, 

1 

but would be difficult to use for a practical system. 

B. Direct Metabolic Rate Monitoring 

If the metabolic rate could be directly determined, it would 

provide an excellent Input for an automatic controller. Once the 

characteristics of a particular LCG have been determined, the 

quantity of heat which must be removed by the LCG to maintain thermal 

equilibrium at any metabolic rate would be precisely known and the 

0 

time constant for response is quite small (approximately 30 s) . 
However, the precise determination of metabolic rate using equipment 
suitable for integration into a Portable Life Support System (PLSS) 
is very difficult. 

A controller using the oxygen consumption rate as a measure of 

0 

metabolism was investigated and the performance was satisfactory, 
but the equipment required to implement this approach for a practical 
system is still unavailable. This is also true of systems that 
measure the CO 2 production rate. 

The rate of decay of pressure in the PLSS 0^ supply could also 

4 ^ 

be used to indicate metabolic rate. However, this is also 
considered impractical since the rate of decay is not very sensitive 
to relatively small changes in metabolic rate and factors such as 
suit leakage cannot be easily accounted for. 

L. Sweat Rate Monitoring 

If the instantaneous sweat rate could be reliably determined, it 
would provide an excellent control signal. The sweat rate responds 
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quickly to changes in metabolic rate and it provides a d"" '•♦'-t 

indication of the subject's thermal state. That is, ’ a s« rates 

are a clear indication of a need for cooling. A swe^' . rate 

(including respiratory evaporation) of about 50-100 ml/h has beer. 

4 8 

suggested as appropriate for thermal comfort, ’ while sv-^ir rates as 
high as 4 1/h have been reported^ during severe thermal stress. 

There are two basic ways of determining the sweat rate of a 
subject who is exercising in an Impermeable suit. The local sweat 
rate from small skin segments can be measured and used with 
appropriate weighting functions to calculate the whole body sweat 
rate, or the total moisture evaporated into the ventilating air 
stream caun be determined by measuring the air flow rate and the 
change in specific humidity. 

The first method suffers from the requirement for 

instrumentation of the subject and the lack of reliable humidity 

cells for use in an operational system. The second method, as 

2 

explored by Chambers, would be well-suited for an operational system 
if a reliable, compact humidity measurement system was available. 
Currently available systems are not sufficiently developed for use in 

Q 

an operational EVA system. 

D. LCG Heat Transfer Rate Monitoring 

The actual quantity of heat removed bv the LCG can be easily 
determined by measuring the water flow rate and the difference 
between the inlet and outlet temperature (AT). Since the 
relationship between heat removal at equilibrium and the total 
metabolic rate can be determined experimentally, the LCG heat 
transfer rate can be used to indicate the total metabolic rate. 
Similarly, the time rate of change (time derivatives) of the LCG heat 
transfer indicates the future trends in metabolic rate. Since the 
water flow rate will be constant in an actual EVA system, information 
on the LCG AT can be used as a controller input. Reliable 
temperature transducers are readily available and no subject 
instrument, 'ion is required to implement such a control system. 
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Control stability can be further Improved if an additional 
temperature measurement, such as heat sink outlet temperature, or 
some subject temperature is added. 

5 

The LCG AT concept investigated by Kuznetz has been selected 
for further development since it requires no subject instrumentation 
and is therefore well-suited for use in aui opei itlonal EVA system. 

Although the instrumentation required to provide the controller 

inputs is very simple, a serious obstacle related to the long time 

constant of the LCG/subject system must be overcome. The delay 

between changes in metabolic heat production and the establishment of 

4 

a new equilibrium LCG AT can be as much as 1 h, ’ this requires the 
use of a sophisticated controller transfer function that 
essentially predicts the future LCG heat transfer needed for comfort 
based on trends in the current data. 

This control scheme has been tested by Kuznetz both 
analytically and experimentally using human subjects. Howev“- the 
human tests, which were intended only to demonstrate syst 
feasibility, were done in an open loop, semi -manual mode. That is 
the instantaneous controller input temperatures were determined and 
manually entered into a computer which then calculated the new LCG 
inlet temperature set point and displayed it to the experimenter. A 
manual valve was then adjusted to obtain the desired water 
temperature. 

A variety of metabolic rate profiles, including some based on 
actual EVA experience, were then performed. The experimental results 
showed that this control scheme was able to maintain the subjects in 
thermal equilibrium for the full range of expected metabolic rates. 
Body heat storage and sweat rates were commensurate with comfort, and 
the subjective reports also indicated thermal comfort. 

The program discussed in this report combined the results of 

5 

this exploratory work by Kuznetz with modern adaptive control 
techniques that can be implemented on a microprocessor. A closed 
loop LCG controller was successfully developed and tested. The 
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present program uses the water teDq>eratures previously discussed as 
the input parameters. However, the adaptive control logic that has 
been developed could be altered to use inputs from humidity or CO^ 
sensors if suitable transducers are developed in the future. 

Adaptive techniques to estimate the transfer function between 
the LCG inlet teiq>erature, and the difference, AT, between and 
outlet temperature, T^, were investigated. This transfer function 
can change with time; therefore, the controller must be able to track 
the changes and make appropriate adjustments to the control based on 
the changed trauisfer function. A brief technical description of 
Kuznetz's previous work in this area and how it relates to the 
present project follows. 

Figure 1 is reproduced from Kuznetz. The data that are plotted 

were obtained from en^irical measurements and from a highly regarded 
•3 

computer model that is in use at the NASA Johnson Space Center 
(JSC). The horizontal axis gives the metabolic rate in watts. The 
left-hand vertical axis gives the LCG heat removal rate in watts when 
1C9 1/h are flowing through the LCG; the right-hand vertical axis 
giv AT. The five curves eU"e plots of AT as a function of the 
metc-oolic rate for a constant T^^. The numbers near points on the 
curves indicate the stored body heat in kilojoules (Hatthours) for 
the indicated metabolic rates. 'iTiere is a well-established comfort 
zone that relates a comfortable body heat storage with a givsn 
metabolic rate. This comfort zone is superimposed upon the five 
curves. All curves are for steady-state conditions. 

It is possible to go frvom a AT on the right-hand axis across to 
a point on the curve at the center of the comfort zone to determine 
an inlet temperature that will keep the subject within the zone. The 
that is necessary to stay at the center of the comf ;rt zone can be 
plotted against AT, as shown in Figure 2 (also reproduced from 
Kuznetz). This curve is very nearly a straight line. 

During steady-state conditions, AT will be pivjportional to the 
metabolic rate because the heat that is removed is proportional to 
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FIGURE 1 PREDICTION OF LCG HEAT REMOVAL (or LCG AT) AS A FUNCTION OF 
METABOLIC RATE. The numbers adjacent to the lines of constant LCG inlet 
temperature represent total body heat storage predictions in kilojoules (Watt hours) 
at metabolic rates of 146, 234, 352, 469, and 686 W. Comfort zone** body heat 
storage vs metabolic rate. 
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the difference between the inlet and outlet water temperatures. If 
there were no dynamics involved, could be adjusted for a given AT 
(or a given metabolic rate) to fall on the curve. In reality, the 
dynamics are highly complex so that a simple change in produces 
complex variations of AT even at a constant metabolic rate. Computer 
experiments with the NASA JSC model were performed to begin to 
quantify these effects. When T^ was moved to fall on the curve of 
Figure 2, the AT moved away from the curve with a new AT. When 
successive new values of T^ were chosen to move back on the curve, 
the system appeared unstable. 

By similar procedures Kuznetz has attempted to derive the change 
in AT versus the size of a step change in T^^. This allows the 
construction of a first-order model. He used the same technique to 
derive coefficients for a smaller second-order model. A second-order 
differential equation was programmed and tested with a subject in a 
LCG. In these experiments, T^ and AT were measured. A new T^ was 
entered into the differential equation and a new AT was obtained. 

This process was iterated until a AT and a T^^ were obtained that fell 
on the curve in Figure 2. The T^ that was indicated to be correct by 
the solution to the differential equation was set by hand into the 
LCG. 

We also used a modeling technique in the SRI study; however, our 
technique is more general and more rigorous and will converge to the 
comfort zone more quickly. In addition, we developed a system to 
control a differential valve that shunts water past the heat 
excheuiger (to control the inlet temperature) with the con 5 )uter for a 
fully automatic ten^erature control. 

In this project, we use well-known adaptive control techniques. 
Figure 3 is a block diagram that illustrates the concept. We assume 
that the LCG (the garment and the enclosed subject) can be modeled as 
a linear system, and that the model can be implemented as a computer 
subroutine. The model input is T^^ and the model estimates AT based 
only on the past values of T^. A subroutine or algorithm adjusts the 
model's parameters to minimize the average of the square of the 
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FIGURE 3 BLOCK DIAGRAM SHOWING THE COMPUTER MODEL 
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difference between the estimated AT and the actual AT. If the model 
were perfect, a test T^ could be entered into the model by the 
computer and the model could be run ahead in time to determine what 
the AT would be for that T^. If the combination did not fall into 
the comfort zone, other values of T^ could be chosen until a 
combination was found that did fall into the zone. The LCG inlet 
temperature would then be adjusted to equal the successful test T^^. 

In practice, of course, it is not possible to make a perfect 
model, but we show in later sections that it is possible to implement 
one that has sufficient accuracy. The iterations to pick a new T^^ 
can be done in several seconds so that this type of control can be 
done in "real time." 

The following were the objectives for a first year of the study: 

• The hardware for a computer driven LCG temperature control 
was to be designed and constructed. 

• There was to be the preliminary development of a computer 
algorithm to control the LCG temperature. 

• There was to be limited testing on human subjects. 

These objectives have oeen met, and the developments will be 
described in the remaining sections. Section II contains a 
description of the experimental system including the cou^uter 
control. Section III contains an overview of the computer program 
that will explain the program's orgamization and design philosophy. 
Sections IV and V explain in detail the valve control algorithm, the 
adaptive model, and the algorithm that determines T^^ with reference 
to the model. Sections VI-VIII describe the experimental results, 
reconmendations for further work, and the Appendix contains a listing 
of the program. 
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II EXPERIMENTAL SYSTEM DESCRIPTION 

Figure 4 is a schematic of the experimental LCXJ temperature 
control system. There is a constant flow through the LCG (nominally 
109 1/h), and water is shunted around the heat exchanger by the 
differential valve; if T^ is to be lowered, less water is shunted. 
Three water temperatures (T^, T^, and the heat exchanger output, T^) 
are A/D converted and entered into the computer. 

An output from the computer is D/A converted to form the analog 
input to an analog servocontrol for the valve position. The three 
temperatures are processed by the computer to calculate a "command 

inlet, temperature", T , which the computer has determined to be the 

c 

appropriate temperature for subject comfort. The computer adjusts T^ 
toward T^ with a control algorithm that includes a closed loop 
control through the computer. (The analog valve servo is in this 
loop.) 

The computer is a Radio Shack TRS-80 Model I. It was chosen for 
the following reasons. We wished to demonstrate the concept with the 
microprocessor-based computer (preferrably a Z-80 microprocessor) to 
show that a microprocessor has sufficient speed (a flight system 
would use a microprocessor) , and that memory requirements would be 
reasonable. Commercial multichannel A/D and D/A converters that plug 
into the TRS-80 are available, and there are commercial temperature 
probes that are compatable with the A/D converter. 

A BASIC compiler for the TRS-80 Model I which compiles programs 
that also run in the interpretive mode is available from Microsoft. 
Debugging is greatly facilitated wtien the program can be run in the 
interpretive mode prior to compiling. The compiled program runs an 
estimated five to ten times faster than the interpreter program. Our 
direct experience when comparing the two versions of the same program 
showed that the interpreter version is too slow by a factor of four. 
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A. The Mechanioal Subsystem 

The breadboard system schematic is shovm on Figure 5. The 
laboratory pumping units are each capable of supplying up to 
3.8 1/min at amy ten^jerature from 5°C to 40 ®C. The heat exchangers 
are brass, tube-in-shell, commercial units puchased from Exergy 
Incorporated, (Part No. 012-122). Standard 3/8 in aluminum tubing. 
Swage fittings, and connectors were used for all plumbing 
connections. A variable speed Micropump was used to set the LOG flow 
rate at 1.9 1/min. (240 Ib/hr). Rotauneters were used to measure the 
flow rate. The flow diverter valve is the saune one used by Chambers 
(Ref. 2). 

This system, including both heat exchangers, was used for the 
initial checkout testa described elsewhere. It was then modified for 
use with an actual LCG, as shown by the dashed lines. The LCG 
replaces the right heat exchanger and pumping utiit. 

The initial human tests were performed using the left pumping 
unit and heat exchanger as shown. However, it was found that, due to 
the approximately 5°C temperature drop across the heat exchanger, it 
was not pc-sible to reduce the minimum LCG inlet temperature to the 
desired 5 to 7°C. This heat exchanger was then bypassed as shown by 
the dotted lines on the left side of the schematic. The circulation 
pump in the laboratory pumping unit was disconnected and the 
Micropump was used to pump through the system. 

In general this system worked satisfactorily. The primary 
problems encountered were leakage and slow response of the diverter 
valve. The effect of this on the test results is discussed 
elsewhere. The large mass of water in the pumping unit's reservoir, 
about 3 1, must change temperature when the diverter valve position 
is chcuiged. This introduces a thermal time constant into dynamic 
performance of the system that is quite different from an actual 
PLSS. It would be highly desireable to go back to the original 
system configuration with the addition of a more effective, 
aerospace-type compact heat exchanger for any further work. This 
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FIGURE 5 BREADBOARD SYSTEM SCHEMATIC 






would more accurately simulate the transient thermal response of an 
actual PLSS heat sink. 

B. The Control Electronics 

Figure 6 is a block dieigram of the control system electronics. 
The water temperatures are measured by silicon diode temperature 
probes that are Inserted Into the water circuit. A l6-channel, 8-bit 
A/D converter is connected to the TRS-80 to sample T^, T^, T^ and AT. 
A measurement resolution of 0.2°C is sufficient for the first three 
temperatures, however, the resolution of AT should be better than 
0.05°C. 

Eight-bit conversion allows a resolution of 0.l6^C over a 40^C 

range for T., T and T. . To increase the resolution for AT, the 
1 o h ’ 

difference between the analog outputs of the T^ and T^ probes, is 
developed with an operational amplifier. The amplifier output is A/D 
converted to 8 bits. The resolution is 0.04°C over a 10°C range. 

An analog servo controls the position of the differential valve. 
The input to the ainalog servo is from a 12-bit D/A converter, and the 
input to the converter is the valve position command from the 
computer. Figure 7 is the schematic of the analog control. Valve 
shaft position as sensed by a potentiometer is fed back through an 
active equalizing network. The 12-bit digital word is output from 
the computer as a sequence of three 4-bit words (most significant 
bits first) that are sequentially stored in a 12-bit latch (see 
Figure 8), prior to the D/A conversion. The servo system can drive 
the valve over the complete range in 5 s. 
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FIGURE 6 ELECTRONIC SYSTEM BLOCK DIAGRAM 
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FIGURE 7 ANALOG SERVO SCHEMATIC 
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Ill COMPUTER PROGRA iRVIEH 


Figure 9 is a sioplified flow diagram of the system cos^>uter 
program. The following major functions are performed: 

• Ten?)erature3 Tj^, T^, and T^^ are measured once a second. 

• Mew valve positic^is aj?e calculated (usually once every 
second) and the position commands are transmitted to the 
analog servo. 

• Model coefficients are up-dated (usually there is a model 
adaption every 20 s). 

• When the operation is out of the comfort zone new T values 
are picked to move back into the zona. 

After this overview, separate sections will be devoted to these 
last three functions. 

Refer to Figure 9 and assume that the program is in the "has 
clock changed" block at the top of the diagram. The program will 
loop back through the block every few milliseconds until the program 
clock changes (it changes once a second). 

At the change the program proceeds down the chart and first adds a 
random step or "dither" to the commanded inlet temperature This 
dither, which will be discussed in detail later, is necessary for 
correct model adaption. One degree C is added to or subtracted from T^; 
whether it is added or subtracted is changed randomly every minute. The 
LCG subject usually will not notice the resulting change in 

Next the three temperatures, T^, T^, and T^, are read. The 
program loops through this box once a second, which is the period 
between rea^-.ngs. A known cuid constant period between samples is 
required for the control of T^^ and for the adaptve model. 

Next the valve position is calculated based on the current 
values of T^, T^, T^, and T^. During this calculation the fraction, 

F, of the water that is shunted around the heat exchanger is 
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calculated. The motor-driven valve requires approximately 5 3 to 
travel from full-shunt to no-shunt, so after a new is calculated 
by the T calculation subroutine, F is held constant for 5s. Then 
feedback through the computer is used to adjust F to drive T^ as 
close as pssible to T^. The F for different valve positions has been 
measured and a non-linear regression curve that is fitted to the 
experimental data is used to determine the valve position. The 
computer ct^smiands this position through the interface. 

There are two possible CRT displays that can be selected from the 
computer keyboard. One is a numerical display of the current values of 
T^, T^, T^, and AT, and the other is a graphics display of (AT, T^) 
and the comfort zone. If the first display is in use it is updated 
every 1 s immediately after the val've command has been given. 

There are four major functions that are performed by the program 
that is depicted by the remainder of the flow chart. These functions 
are: model estimate of AT, model adaption, CRT display update, and 

T^ calculation. They are done once each in an interval that usually 
is 20 3 long (the length is set by the parameter L0%). Each time the 
program loops through the counter J0% is incremented until J0% = L0%. 
Then JOJ is set to 0 (modulo LOU). 

The first two functions are done during the same 1 s long 
interval, and the display and T^ calculations are done during 
separate subsequent 1 s intervals. Because the total computation 
time for all four functions is greater than 1 s the program would not 
always be able to loop through in 1 s if all the calculations were 
done together. 

When J0% ^ 1. a new filter (model) output, AT, is calculated. 

Thia need be done no more often than every 20 s because empirical 
measurements have shown that this period is sufficient for the 
bandwidth of an actual LCG (the sampling rate shoulo be greater than 
twice the actual bandwidth). The experiments that determined thac 
L0% = 20 is appropriate are described later. A new update in the 
estimate of AT is necessary for an update in the model adaption, so 
the adaption must follow the estimate of AT. Both functions will be 
described in detail in a later section on modeling. 



When J0% = 2 the graphical display is updated if it is in use. 
This display is a real-time equivalent of Figure 2 where the latest 
"target" (AT, T^) is plotted along with the latest actual (AT, T^^). 
These two points are plotted along with the boundaries of the comfort 
zone, so that system performance can easily be followed. 

When J0% = 4 the most complex and time consuming of the 
con^jutations begins. First there is a test to determine if there 
have been a sufficient number of 20 s or LOJl intervals since the last 
calculation. The variable J4J is incremented by one; unless it is 
greater than or equal to parameter SH% the program branches directly 
back to the main loop so that there is no new calculation of T^. 

There must be this delay in the calculation of the next T^ 

because the system cannot move into the comfort zone immediately. If 

a new T^ is chosen before the system has reached steady-state within 

the zone, the calculated T^ will be pushed farther from the correct 

c 

T by each new calculation (i.e., the system will be unstable). The 
c 

parameter SH$ is chosen so that the "dead-time" is 4-5 min (SH% is 
12-15). This range atnd other techniques that will be discussed 
result in stable operation. 

Next it is determined whether or not (AT, T^^) is outside of the 
t'orfort zone. If not, the program returns to the top of the flow 
diagram to await the beginning of the next 1 s interval. If it is 
outside, the T^ calculation is started. This calculation may last 
more than 1 s; if it does a flag indicating that the calculation is 
not complete is set and the program returns to the top of the loop. 
When the program passes through the "does J0% = 4 test" (J0% now is 
greater than 4 ) , a flag-set indication will cause a branch to the 
calculation portion of the subroutine and the T^ calculation will 
continue. Thus, there will still be a pass-through the main loop 
once a second. When the calculation is con^lete the flag is reset so 
there is no branching on subsequent passes during the 20-s interval. 
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Figure 10 gives more of the calculation details. On entry to 

the calculation portion the "is flag set" is tested. The answer "no" 

means that the calculation has not started, so an initial trial value 

for T is chosen. After a teat to check that the 1 a interval has 
c 

not ended, a model steady-state AT (the AT that the model predicts 

will occur in 5-6 min if T, is set to the value for T ) is 

1 c 

calculated. 

A test is made to see if the model (AT, T ) xs in the comfort 

c 

zone* If it is, there is a branch through the reset box and back to 

the main loop (return)* If not a counter (J8%) is incremented to 

indicate how many trial T^s have been tested. If the number is less 

than L8% (usiaally set to 10 or 12) the program returns to the "next 

trial T " box* 
c 

Different procedures for picking the "initial" and "next" values 
are used when the main loop with the flag set the "next trial T^" 
procedure is used* 

If a is not chosen that results in (AT, T^) being in the 

comfort zone by the (L3% + 1) pass, a single pass calculation 
technique is used* In general, the actual AT will not agree with the 
predicted AT as closely when this latter technique is used; however, 
a solution is guaranteed. After the calculation, the program returns 
to the main loop* Both the multi- and single-pass techniques 
calculate AT using the adaptive model and the current values of the 
coefficients* 

The next three sections will describe in more detail several of 
the algorithms that are used in the subroutines* 
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FIGURE 10 SIMPLIFIED FLOW DIAGRAM OF Tq CALCULATION 
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IV VALVE CALIBRATION AND SERVOCONTROL 

A simple valve calibration technique has been developed. Since 
the theory of the calibration is related to the servocontrol theory, 
both are described in this section. 

A. Valve Calibration 

Assume that there are W grams of vater at a temperature of Tj^°C. 
When a fraction, F, of this water is heated to T^°C the number of 
calories that is added to the fraction is 


AC = F*W*(T - T„) 
o n 


( 1 ) 


Now mix the water that has been heated to T with the remainder of 

0 

the water (that still is at a temperature T^) . The rise in 
temperature, AT, of the mixture is 


AT 


W 


( 2 ) 


In the LCG system 


T. r at + T. 


F-W-CT -T^) T^ 

= - Q- + h 

W 

= F(T^ - T. ) + T. 

on n 

or 


F = 


r - T 

h jl 

T - T 
o h 


( 3 ) 


( 4 ) 


Thus it is possible to measure the fraction of water that is 
shunted around the heat source by measuring and T^. During 
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the calibration the valve position was stepped from full shunt to 
zero shunt and F was calculated from Eq. 4. A non-linear regression 
curve of valve position as a function of F was fitted to the data. 

The curve and data are shown in Figure 11. The curve is nearly 
linear and not correct when F is almost 1. This valve did not allow 
a total shunt; there was always some water going through the heat 
exchanger due to internal leakage. This annoying feature degraded 
the still satisfactory servo performance. 

The equation of the curve is given in lines 7090 zuid 7110 in the 
listing (Appendix) where CP(NJ) are the coefficients. The values are 
given in line 35000. 

B. ServoControl of Inlet Temperature 

We have designated T^ to be the inlet temperature that is 
designated by the computer to reach the comfort zone. Let F' be the 
shunt fraction that gives an inlet temperature of T^. Let F be the 
present fraction that resulted in the present temperature, T^. Then 
using Eq. 3 

T, = f ' < W‘V * "h * • <5) 

where the new T. is (T. + AT. ) , 

h h h 

- II . (F'(T^ - - A-n_) 

+ t + AT.) 

n n 

-[F'(T^ - T^) 

= (F' - F)(T^ - T^) 

+ (I - F') ATj^ (6) 
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where AT^ is the change in with F ' . 
Neglecting 


F' 



( 7 ) 


Equation 7 is used to set the initial value of F' whenever there is a 
change in T^. This fraction will not be exactly correct so feedback 
through the computer is used for further adjustment of F* after the 
initial value is set. 

Approximately 5 s (a 5 s dead time) are allowed for the motor to 

drive the valve to the correct position before feedback through i.he 

con?)uter is used to account for AT. i 0 and other errors. A first 

h 

order difference-equation equalization filter was programmed and 
placed in this loop in an effort to speed the response of the 
through-the-computer loop. This equalization filter was abandoned 
because of an increase in loop noise; the difference equation is an 
imperfect differentiation. 

Instead of using equalization filters the following rules are 

used: 

If T - T, > 0 

c i 

F(n) = F(n-l) + C^(T^ - T^) (8) 

If T - T >1 
i c 

F(n) = F(n-l) - (9) 

If 1 > T - T >0 
i c 

F(n) = F(n-l) + € 2 ( 1 ^ - T^) , (10) 

where F(n) is the fraction at the n‘ ‘ second after the end of the 
dead time. 
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Note is the current value and each second a new F(n) is calculated 
using one of the three rules. Equation 3 is used when the command 
temperature is greater than the present inlet temperature. The 
increase in F is proportional to the difference. 

Equation 9 is used if is less than the present value of 
and the difference is greater than I^C. Equation 10 is used when 
must go lower but the difference is less than l^C. Equations 9 and 
10 allow a maximum change in F of C2 each second; but when comes 
near to the change is proportional to the difference. 

The inlet temperature can be changed rapidly in the negative 
direction because the water is quickly cooled by tne heat-exchanger. 
If Eq. 10 is used over tne whole range when T^ is regative-going 
(i.e., Eq. 9 is not used) either the change in F will be too large 
and there will be unacceptable overshoots or must be made small to 
damp out the overshoots. If a sufficiently small is used to damp 
the overshoot the convergence will be too slow when T^ is close to 

T . 

c 

Because the only source of heat is from the subject, T^ can be 
lowered much more quickly than it can be increased. This is why 
there are two gain constants, and C^. Their values were 
empirically determined to be = 0.075 and = 0.02. They depend 
on many system parameters such as flow rate, water volume, and LCG 
heat transfer characteristics. Because of all the unknown or hard- 
to-measure parameters it probably is best to determine and by 
experiment instead of by engineering design. 

The above algorithm may appear to be a "seat-of-the-pants” 
approach and to a certain extent it is. Unfortunately, classical 
servo design techniques are far from applicable because of system 
non-linearities and most importantly because there is no source of 
hot If the servocontrol could rapidly dump a large number of 

calories into the water loop, then T^ could rapidly be driven in both 
positive ana negative directions. It then would be a nearly 
classical servo design problem, and equalizing filters (difference 
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equation iterated by the computer) to maximize the rise time and 
minimize the oversnoot and settling time could be used. Optimum 
control techniques could be used to minimize the time required to 
return to the comfort zone. 

After many months of intense experimentation the algorithm has 
proven to be satisfactory, and given the system constraints 
(especially the valve characteristics) we believe that there is 
little likelihood that the design could be significantly improved. 
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V THE MODEL 


A. Model Deaorlption 

The model (implemented in a subroutine) Inputs are measured 
values of T^. Not all T^^ are used; only every L05tth (usually every 
20th) T^ is input. It was mentioned earlier that this low sampling 
rate can be used because of the low LOG bandwidth. The model outputs 
are the predicted values of AT, AT* at the times there is an input 
(i.e., every L0$ seconds). Model parameters are adjusted (adapted) 
in real-time to give the minimum mean-squared difference between AT 
and AT*. 

The model has been tested against a number of subjects in an LOG 
and against data from the 41 -node LCG model at the JSC. The model 
generated AT* closely follows the LCG and the 41 node model 
calculated AT. 

The 4 1 -node model is a complex simulation. The parameters in 
our less complex model do not represent different parts of the 
anatomy or heat flow from the anatomy to the LCG. Thus our adaptive 
model cannot oe used to directly indicate body heat stored or heat 
flow from a portion of the thigh. However, it requires little memory 
and calculacion and can be implemented on a microprocessor to run 
much faster than real-time to predict the important variable, AT. 

The model alone is a definit ' contribution of the project that may 
have other applications. 

Figure 12 is the model block-diagram. The very-low-frequency 
components of T^^ are passed through a low-pass filter, multiplied by 
a gain of (C^-1) and added to the adaptive filter output. The high- 
frequency components of T^^ form the inputs to the adaptive filter. 

The logic of this high-pass, low-pass scheme is as follows. The 
waveform of the temperature as a function of time has a large, slowly 
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FIGURE 12 MODEL BLOCK DIAGRAM 
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varying average (or oC) omponent and a muon smaller and more rapidly 
varying high-frequency component. During the examination of 
experimental data it was noted that the average of was almost (but 
not exactly) proportional to the average of AT. The DC value of AT*, 
is inserted through ' he low-pass gain path. Then the adaptive filter 
has to adapt only to the much smaller AC component, and the adaption 
accuracy can be leas without degrading the overall performance. 

Both the low-parc and high-pass filters are realized with first- 
order difference equations.* The low-pass equation is (see line 6105 
in the listing): 

V"> = Vl<“> * • <"> 

Where A seta the bandwidth of the filter (usually A s .05), T.(n) 
m iii 1 

is the nth input T^ amd 1^(0) is the filter output. The high-pass 
filter is 

H (n) = T,(n)-T,(n-1)+(1-B„)H.(n-1) , (12) 

1 11 mi 

where Hj^(n) is the filter output and B^ sets the low-frequency out 
off (usually B^ = 0.01). Both A^ and B^^^ were detern.lned empirioaliy. 

The variable C^ (see Figure 6) is approximately equal to the 

ratio of the average of T to the average of T. . The variable I (n) 

o 1 m 

Is the average (mean) of T^^(n). An equation of the same form as Eq. 

11 [except that T^(n) is substituted for T^(n)] is used to estimate 

the mean of T , 0 (n) (line 6105). Since 
o m 


the mean of AT, D^, is 


*A sometimes awkward variable notation is used in the equations in 
this report. The notation will be :imilar to the variables used in 
the computer program to help with the program documentation. 
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D„(n) s 0(n) - I (n) 

In ID ID 

» (C.-1) . (13) 

IB <1 

A formula for has baen found axparlmantally that is more 
accurate than the ratio of I„(n) to 0_(n). It Is (line 9115) 

ID ID 

0 (“) r -} 

'd* rwK‘V">-vJ • <"■> 

m 

where is a parameter that is usually set to 0.017. The equation 

was obtained by fitting a regression curve through points of a plot 

of C. verses (I„(n)-T„). 
d me 

Figure 13 is a schematic of the adaptive filter. This form is 
sometimes called a tapped delay line or finite impulse response 
filter (FIR). In the computer the output is given by 

m 

D (n) - 2 H (n-J+I) C . (15) 

J-l ^ 

where the Cj are the adapted coefficients. 

A discussion of the theory of an FIR is beyond the scope of this 
report. However, if m is made sufficiently large the FIR can be made 
to simulate any linear system with one input and one output channel. 

Also, there exist well-understood algorithms to adjust the to 
minimize the mean-squared difference between AT and AT*. Ours is a 
slightly fflcdifiod version of the Wldrow-Hoff algorithm. In our 
algorithm 

Cj(n+1) a Cj(n) ♦ kgA^H^(n-J+1) , (16) 

where Cj(n) Is the value of at the nth sampling time. 

If (AT- T*)^ > 1, then 
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FIGURE 13 SCHEMATIC OF THE 
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sign of (AT-AT*) 


(17) 


If (AT-AT*)*^ s 1, then 


k = (AT-aT*) , (18) 

c 

where A| is an adjustable constant; A^^O.OI. 

During each step the change in the weight at the Jth tap is 
proportional to the value cf at that tap times the error between 
AT and AT* if the error magnitude is <»ie or less; it is proportional 
to ±1 if the magnitude is greater than one.* Notice that only two 
input variables are required: AT and 

Examination of Eq. 16 shows that if T^ is constant (and 
therefore = 0) there will be no adaption. To ensure that adaption 
is occurring continuously the small random step or dither is added to 
T^ 30 that Tj^ is never constant. The size and frequency of change of 
the dither depend on the LCG bandwidth euid temperature and other 
system noise. A step size of I^^C has proven to be more than 
adequate (probably a smaller size will do) , and choosing a random 
value of the square-wave once a minute is adequate. 


B. Model Experiments 

Upon completion of the control program, a series of human 
subject tests were conducted to refine the LCG predictive model and 
adaptive subroutine coefficients. As discussed earlier the purpose 
of these subroutines is to provide a continuously updated predictive 
model of the LCG transfer function. This allows the LCG performance 
to be predicted before the valve position is actually changed. Thus, 
the desired system response can be selected based on the current 
state of the system. 


*The truncated error function allows stability with a much larger 
Af and therefore a much faster adaptation. 
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A set of five model coefficients must be used to minimize the 
error between the actual LCG AT values and these predicted by the 
computer control algorithm during a given exercise condition. This 
set of model coefficients was determined in the following manner: 

• The control system was operated without T(. calculation during 
treadmill exercise tests with specified LOG inlet 
temperatures. The cooling system's transient response 
characteristics were obtained by recording the resultant LCG 
outlet teiq)eratures during these tests. 

• A parametric computer study was then performed to determine 
the set of model coefficients that could best be used to 
euialytically duplicate the above experimental results. 

If the calculated LCG outlet temperatures were sufficiently like 

those that were obtained from the actual tests, the model 

coefficients could be used to calculate actual LCG outlet 

temperatures ahead in time with sufficient accuracy for use in the 

automatic controller. 

Four subjects were tested in this series of experiments. In 
each test sequence, the LCG flow rate was set at 1.9 1/min 
(240 Ib/h). During each test the subject wore insulated clothing, 
including gloves and hat, over the LCG while he performed the 
experiment. Bach test sequence consisted of a 10-min standing rest 
period, a walking treadmill exercise period of approximately 20 min, 
and a standing recovery period of 10 min. An inltl 2 d treadmill speed 
of 3 mph (metabolic requirement, = 800-900 BTO/h) was used 
during each exercise period. This speed was held constant during 
three of the tests. During the fourth test the treadmill speed was 
later increased to 5 mph * 2400 BTU/h). 

The system shown on the schematic in Figure 5 was modified for 
these human tests. The actual LCG was used in place of the heat 
exchanger and laboratory pumping unit heat source on the right. It 
was also found that due to the temperature drop across the left, heat 
sink, heat exchanger it was not possible to reduce the minimum LCG 
inlet temperature to the desired 40-46°F range. The system was 
modified so that the flow diverter valve would send fluid directly 
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through the pumping unit rather than through the heat exchanger. 

This allowed the heat sink outlet temperature to be brought down to 
the desired range. The conq)leted system was then ready for the human 
tests. 

The exercise profiles, specified or command LCG inlet 
temperatures, the actual or measured L(X} inlet temperatures and the 
measured LCG ATs for each test sequence are shown in Figures 14, 15, 
16, and 17. The actual inlet temperature, shown in the middle of 
each figure, includes the ±1°C dither that was superimposed on each 
commanded inlet teaqperature. The solid, measured AT trace in the top 
portion of each figure shows the difference between the actual LCG 
inlet and outlet temperatures for each test. 

The inlet temperature profile for each test served as the 
independent Veu*iable during the various pai*ametric calculations. The 
objective of the parametric study was to determine the set of 
coefficients that could best duplicate all of the measured AT traces 
shown in Figures 14-17. 

Table 1 provides the typical results from one of the analytic 
calculations. The various parameters listed at the beginning of the 
printout au*e discussed in Section V of this report. The variation of 
the model as the system adapts to different conditions is shown by 
the changes in the values of the filter coefficients. The actual 
performance of the model is listed as a function of time. The actual 
inlet temperature, which included a ±1®C dither signal in this case, 
is listed in the third column. The resulting actual LCG AT is listed 
along with the filter output, which is the model’s predicted AT. The 
cumulative RMS error between the actual and predicted values is then 
listed after the first 1000 s of experiment time. This error was 
used to help to refine the values of various parameters in the model. 
These parameters were adjusted so as to produce the smallest RMS 
error. The last two columns listed are not pertinent to this 
discussion. 
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FIGURE 14 LIQUID COOLING GARMENT TEMPERATURE CONTROLLER OPEN LOOP 
TEST RESULTS --SUBJECT ONE 
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FIGURE 15 LIQUID COOLING GARMENT TEMPERATURE CONTROLLER OPEN LOOP 
TEST RESULTS --SUBJECT TWO 


40 









SPEED - MPH TEMPERATURE - °C TEMPERATURE 




ELAPSED TIME - sec 


FIGURE 16 LIQUID COOLING GARMENT TEMPERATURE CONTROLLER OPEN LOOP 
TEST RESULTS --SUBJECT THREE 


41 



SPEED - MPH TEMPERATURE - °C TEMPERATURE 



ELAPSED TIME - sec 

FIGURE 17 LIQUID COOLING GARMENT TEMPERATURE CONTROLLER OP«=N LOOP 
TEST RESULTS--SUBJECT FOUR 
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Table 1 

SAMPLE ALGORITHM PROGRAM OUTPUT 


OF.ivj:;N*M PAGE ?S 

OF POOR QUALITY 


THE HUriBEP OF ThPS IS 4 , THE SIZE OF THE DITHEP IS 2 DEG. C. 
-JTWE. -NUMBER OF TEN SECOND INTERNALS BETNEEN CALCULATIONS IS 2 . 
THE ADAPTION IS NOT HELD. THE FILE NAME IS DAUE. 


THE ADAPTIVE CONSTANT IS NOW .05 . N IS NOW 0 


THE VALUE OF THE MEAN CONSTANT NOW IS .2 . N IS NOW 0 


N® 3^ . THE FILTER COEFFICENTS ARE AS FOLLOWS: 
-.378675 .11631? .115A68 .1287^1 


N= 64 . THE FILTER COEFFICENTS ARE AS FOLLOWS: 
-.64267B -.18342 .121965 .172501 


N*_24— . -IHE-EILXER COEFFXCENTS ARE AS. FOLLOM21 

-.838757 .0988044 .119488 .168583 


N 

TIME 

INLET 

DELTA 

FILTER 

RMS 

T. 


DC 


(SEC.) 

TEMP. 

TEMP. 

OUTPUT 

ERROR 




5 

190 

20e4222 

2.7027 

-'13.6021 

0 

16 

12 

1.20935 

6 

120 

20.7735 

1.99089 

3.13627 

0 

16 

32 

1.18011 

7 - 

-14-0 

- 19.594 

2-* 6923Z 

- -J1.931X5L- 

0_. . 

16 

52 

1.17016 

8 

160 

19.7835 

1.75467 

2.83041 

0 

17 

12 

1.15165 

9 

180 

Z0.1232 

1.62464 

2.53514 

0 

17 

32 

1.13576 


200 

-ZO^liU. 

1 .6691 3 

- -2.28304 - 

n 

17. 


1 . 1 

11 

220 

18.5924 

3.03786 

2.05846 

0 

18 

12 

1.13213 

12 

240 

18.7481 

2.72648 

2.2184 

0 

18 

32 

1.13483 

13 

260 

19.0994 

2^ ai46B 

7 1 79^2 

0 

18 

52 

1.12861 

1^ 

280 

17.5019 

3.56135 

2.27427 

0 

19 

12 

1.14302 

15 

300 

17.6247 

2.75716 

2.77902 

0 

19 

32 

1.14561 

^1^ 

32M 

_ia. 0534 

2-.449J-9 

- y 

-0 _ 

19 

52 

-1.14364 

17 

340 

17.9093 

2.17881 

2.68415 

0 

20 

12 

1.1393 

18 

360 

18.065 

1.86742 

2.4239 

0 

20 

32 

1.13214 



-380 

18.9669 

1.16252 

2^14418 

0 

20 

52 

1.11745 

20 

400 

19.2667 

1.12152 

1.98435 

0 

21 

12 

1.1051 

21 

420 

19.5397 

1.12494 

1.92642 

0 

21 

32 

1.09517 


-44iJ 

17.6799 

3.20548 

2^28895 

0 

21 

52 

1.11164 

23 

460 

19.111 

1.4329 

2.05032 

0 

22 

12 

1.10411 

2^ 

480 

18.0534 

2.44919 

2.13096 

0 

22 

32 

1.11027 

25 

500 

1748648 

2.26778 

2^00955 

0 

22 

52 

1.11351 

26 

520 

18.0872 

1.82294 

2.02911 

0 

23 

12 

1.111 

27 

540 

18.065 

1.86742 

1.85264 

0 

23 

32 

1.10949 

28 

560 

18.6938 

1.1591 

1^8039 

0 

23 

52 

1.09982 

29 

580 

19.2444 

1.166 

1.54285 

0 

24 

12 

1.09167 

30 

600 

1<5^.2667 

1.12152 

1.64°72 

0 

24 

32 

1.08477 

31 

620 

19.5175 

1.16942 

1.6273 

0 

24 

52 

1.07<^62 

32 

640 

10*5397 

1 .12494 

1.61709 

0 

25 

12 

1.07508 

33 

660 

1*5.795 

1.1729 

1.41534 

0 

25 

32 

1.07181 


43 



ORIGINAL PAGE IS 
OF POOR QUALmr 


N 

TIME 

INLET 

DELTA 

FILTER 

RMS 

T4 

DC 

— 

(S£C^) 

TEMP. 

TEMP . 

OUTPUT 

ERROR 


- - 

3^ 

680 

19.795 

1.1729 

1.40955 

0 

25152 

1.06923 

35 

-70 0 

19.7283 

1^30635 

1-^39763 

0 - _ 

26112 

1.06862 . 

36 

720 

20.0^58 

1.22081 

1.27262 

0 

26132 

1.06703 

37 

7^0 

17.3«^07 

3.78377 

2.32001 

0 

26152 

1.09447 

38 

Tsa 

-J.9.1555 

1.3€3gg- 

l^ji266a _ 

-- Q. 

-271-12- 

L. 

3? 

780 

19.7506 

1.26187 

1.24937 

0 

27132 

1.08431 


800 

19.7728 

1.21739 

1.17135 

0 

27152 

1 . 07966 



820 

19.7061 

1.35084 

-1.53289 . 

0 

28112- 

X.-07741 

^2 

8^0 

19.1217 

1.97019 

1.85307 

0 

28132 

1.08247 

^3 

860 

19.3886 

1.43638 

1.64655 

0 

28152 

1.08079 

£4 

880 

17.^468 

3.11303 

__2. 54576- _ 

--0 

29112 

1.09874 

^5 

900 

18.749 

1.60742 

1.6391 

0 

29132 

1.09617 

^6 

920 

18.6156 

1.87432 

1 . 77264 

0 

29152 

1.09706 

_ _47_ 

9^0 

17.9982 

2.00087 

JL. 96617 

0 

30112 

- 1.09978 

^8 

960 

18.0205 

1 . 95639 

2.03261 

0 

30132 

1.10148 

^9 

980 

18.6716 

1.20358 

1.58136 

0 

30152 

1.09404 

511 

1000 

^8.9669 

_1. 16252.- 

._1a35423__ 

_.1?1713_ 

31112 

—1.08738 

51 

1020 

19.2444 

1.166 

1.28136. 

.158213 

31132 

1.08189 

52 

10^0 

17.9867 

2.58264 

2.15126 

.280568 

31152 

1.09381 

53 

1060 _ - 

_ 18.838 

1.42948- 

-^1.59695 - 

_ ,257003 

32112 

-- 1.09019 

5^ 

1080 

17.6914 

2.62371 

2.29538 

.272765 

32132 

1.10133 

55 

1100 

19.2222 

1.21048 

1.10921 

.252408 

32152 

1.09341 

5A_ 

1120 - 

-.19.2222 

1.21048 

1.46714-_ 

. .25302 

33112 

1.08716 

57 

11^0 

19.5175 

1.16942 

1.21078 

.23713 

33132 

1.08153 

58 

1160 

17.469 

3.06855 

2.74314 

.248492 

33152 

1.0992 

-59- 

1180 

19.0888 

J.. 47739 

1.45888 

.-235813 

34112 

1 . 09475 

60 

1200 

18.1424 

2.27126 

2,18097 

.226481 

34132 

1.10069 

(yl 

1220 

19.2 

1.25497 

1.17564 

.218045 

34152 

1.09344 

62 

i 2 aa - - 

1 Q . ^777 

1^2104a- 

.-1.S678Z- 

..231752 

3511.2 

1.08722 

63 

1260 

19.2444 

1.166 

1.45124 

.235975 

35132 

1.08179 

6 ^ 

128C 

19.795 

1.1729 

1.21775 

.228267 

35152 

1.07712 

65. 

1300 

18.0196 

3.07545 

- -2. 45745 

.269666 

36114 

1.09499 

66 

1320 

17.9199 

2.7161 

2.41207 

.271807 

36134 

1.10584 

67 

13^0 

17.4806 

2.48678 

2.63892 

.266572 

36154 

1.11275 

68 

1360 - - 

16.8632 

--J2.-61333 

2.567Z4__ 

^9678 

37114 

1.12061- 

6 ? 

1380 

16.5457 

2.69887 

2.59695 

.254126 

37134 

1.1285 

70 

1^00 

15.7505 

3.1813 

2.96869 

.252304 

37154 

1,14181 

. 

„ - 1020 

14.8486 

3.8862 

3.4648 

.262366 

38114 

1.1629 

72 

1^40 

14.9598 

3.6637C 

3.49052 

.259129 

38134 

1.17777 

73 

1^60 

14.6822 

3.6603 

3.70517 

.253839 

38154 

1.19078 

7A 

1^80 

13.758 

4.40968 

-4.25Z9 

.230679 

39114 

1.21358 

75 

1500 

13.7696 

3.82791 

4.2967 

.262441 

39134 

1.22518 

76 

1520 

13.5365 

3.73546 

3.93005 

.260244 

39154 

1.23439 



15^0 

13.581 

3.64649 

3.54269 

.256306 

40114 

1.24071 

78 

1 560 

14.1903 

3.51994 

3.14068 

.2615. 

40134 

1.24211 

79 

1500 

16.1916 

1.1725 

1.9029 

.289639 

40154 

1.20522 

HU 

1600 

16.4025 

1.30943 

1.76362 

-^296376 

41114 

1,17815 

81 

1620 

15.8955 

2.33262 

2.43935 

.292318 

41134 

1.17164 

82 

16^0 

15.251 

2.42159 

2.83924 

.296895 

41154 

1,16776 

--83 

1660 

15.78^3 

2.55504 

2..i9597 

.293493 

42114 

1.16<'56 

8^ 

1680 

15.8732 

2.3771 

2.47468 

.28074 

42134 

1.16314 

85 

1700 

16.1907 

2.29156 

2.16911 

.286415 

42154 

1.15869 
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N 

TIME 

INLET 

DELTA 

FILTER 

RMS 

T$ 

DC 


(SEC^) 

TEMP. 

TEMP . 

OUTPUT 

ERROR 



86 

1720 

16.213 

£.2^707 

2.20602 

.282599 

^3:1^ 

1.15456 

- -87 

17^0 

16.2574 

2*15811 

2.2291? 

.279094 


1.1501 

86 

1760 

16.2574 

2*15811 

2.27315 

.276108 

^3:5^ 

1.14657 

89 

1780 

16.2797 

2.11362 

2.23471 

.273306 


1.14317 

9^ 

0804- 

-16.3019 

0694:4 

7.19711 

.270694- 

_ 44:3^ 

1 .13988 

91 

1820 

17.4814 

1.36772 

1.14631 

.269622 

44:59 

1.12671 

Q2 

18^0 

17.5926 

1.1453 

1.22497 

.266746 

45:19 

1.11366 

- 93 ^ _ 

4860 

18.1655 

1610771 

.-836209 

- .-266855 

45:39 

1.10235 


1880 

18.7383 

1.07013 

.743528 

.268327 

45:59 

1.09254 

95 

1900 

19.01t~ 

1.07355 

.880402 

.266918 

46:i9 

1.08476 

-94^-- 

4920 

19.0113 

1^07355 

- 4^24168 

- .2652 

4A:39 

1.07875 

97 

19^0 

18.4047 

1 .73739 

1.77136 

.262468 

46:59 

1.08195 

98 

1960 

19.2222 

1.21048 

1.07038 

"60546 

47:19 

1.07795 

99 

1980 

19.2444 

1.166 

U 13426 

.257967 

47:3? 

1.07433 

100 

2000 

19.2444 

1.166 

1.147_6 

.255438 

47:59 

1.07149 

101 

2020 

19.1777 

1.29945 

1.37478 

.253186 

48:i9 

1.07072 

1412 

-^0^0 - 

-19.1555 

4.34394 

- 1-. 36288 

- .2508 

48:3® 

1.07061- 

103 

2060 

19.1332 

1.38842 

1.36282 

.248491 

• 48:59 

1.07101 

10^ 

2080 

18.1869 

2.18229 

2.10342 

.246451 

49:19 

1.08052 

- -445 

-.2400 

- 18.065 

1.86742 

. -2.089-63 

_ ^246039 

-49:39 

1.08496 

106 

2120 

18.1317 

1.73397 

1.82586 

.244175 

49:59 

1.08705 

107 

21^0 

18.1762 

1.645 

1.57651 

.242228 

50:i9 

1.08773 

- 448 

2160 - 

18.2206 

1.55603 

4-.^5099 

- .240467 

5D 139 

-- 1-. 08727 

109 

2180 

18.7161 

1.11462 

1.18577 

.238334 

50:59 

1.08165 

110 

2200 

19.0113 

1.07355 

1.02055 

.23647 

5i:i9 

1.07649 

Ill 

2224 - 

19.0113 

1.07355 

_ 4.-15344 

-.234775 

51:39 

1.07241 

112 

22^0 

18.4269 

1.69291 

1.73555 

.232966 

51:59 

1.07627 

113 

2260 

18.5826 

1.38152 

1.57915 

.232455 

52:19 

1.07588 

--14^ 

2280 - 

- 17-.9537 

-Z. 08984 

-4^97452 

-- .231403 

_52:39_ 

-2.08376 

115 

2300 

18.7- • 

1.11462 

1.21786 

.229697 

52:59 

1.07887 

116 

2320 

19.03j6 

1.02907 

.995254 

.228014 

53:19 

1.0738 

_ 117 

23^0 

19.3112 

1.03255 

- .840784 

.227924 - 

53:39 

1.06961 

118 

2360 

19.3112 

1.03255 

1.09007 

.226372 

53:59 

1.06631 

119 

2380 

19.5842 

1.03597 

.956094 

.224952 

54:19 

1 . 06355 

-424 

2^40 - - 

49.5842 

1-.-03597 

-- 4.^06875 

- - .223396- 

-54:39 

. -. 2-.06137 

121 

2^20 

17.5135 

2.97958 

2.72701 

.223813 

54:59 

1.08163 
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Table 2 shows the results of variations in 3 of the model 
parameters, with the resulting effects on the cumulative RKS error, 
for four different subjects. The combination chosen based on this 
preliminary optimization was number of Taps s 4, adaptive coefficient 
= 0.05, and mean coefficient s 0.2. This combination produces an 
average error of about 0.4^C for all of the subjects. This error can 
be reduced even further by more detailed refinement using additional 
test data from more subjects. These values were then used to 
generate the various calculated AT results shown as the dashed traces 
in Figures 14-17. In each case, the adaptive coefficients were set 
equal to 0 at the start of the run. This produced a relatively large 
initial error that decreased after 5-10 min as the model adapted. 

This initial adaptation will eventually be avoided when the optimal 
starting values for the coefficients have been identified. The 
figures show that following the initial adaptation, the predicted auid 
actual values are very close, even during large transients. This 
clearly demonstrates the accuracy of this adaptive technique for 
predicting LCG performance. The values of the various program 
parameters derived from this series of calculations were included in 
the automatic LCG inlet temperature control program. 
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Table 2 


RMS - ERROR MINIMIZATION 



ADAPTIVE 

MEAN 


SUBJECT 


// TAPS 

COEFFICIEIiT 

COEFFICIENT 

1 

2 

3 

4 

2 

.15 

.1 

.5119 

.3461 

.4916 

1.366 

2 

.15 

.2 

.4245 

.3041 

.4278 

1.0039 

2 

.05 

.1 

.5180 

.4213 

.4799 

.7220 

2 

.05 

.2 

.4787 

.3991 

.4416 

.6006 

2 

.025 

.1 

.6077 

.4899 

.5834 

.8018 

2 

.025 

.2 

.5864 

.4695 

.5488 

.7596 

4 

.15 

.1 

.6726 

.7606 

.6900 

2.3674 

4 

.15 

.2 

.5657 

.2387 

.3317 

1.6281 

4 

.05 

.1 

.4899 

.2555 

.3466 

.7350 

4 

.05 

.2 

.4619 

.2220 

.3051 

.6454 

4 

.025 

.1 

.5743 

.3418 

.4428 

.7512 

4 

.025 

.2 

.5613 

.3336 

.4206 

.7235 

6 

.15 

.1 

.8443 

1.2502 

1.1426 

2.9721 

6 

.15 

.2 

.7493 

.4733 

.j519 

2.1357 

6 

.05 

.1 

.4417 

.3991 

.4309 

.9090 

6 

.05 

.2 

.4330 

.2251 

.3252 

.6524 

6 

.025 

.1 

.5325 

.3430 

.483« 

.7631 

6 

.025 

.2 

.5332 

.3274 

.4458 

.7248 

8 

.15 

.1 

1.1326 

1.6311 

1.5000 

4.0U5 

8 

.15 

.2 

.8105 

.7234 

.8092 

2.7082 

8 

.05 

.1 

.4728 

.4849 

.4957 

1.0764 

8 

.05 

.2 

.4354 

.2515 

.3378 

.6973 

8 

.025 

.1 

.5453 

.3633 

.4847 

.8242 

8 

.025 

.2 

.5358 

. 3 /") 

.4591 

.7379 
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VI CHOOSING THE INLET TEMPEBATORE 


Aitsuma that the model has adapted and that due to a metabolic 

transient T^) is out of the comfort zone. The computer now 

picks a new T by entering a trial T into the model and running the 
c c 

model ahead in time to see if the model (AT*, T ) falls into the 

o 

zone. If the first choice is incorrect others are chosen iteratively 

until (AT*, T^) falls into tne zone. If there is no success after 

the LSK trial T s, a second more direct (but less accurate) method 
c 

that does not require iteration is used. This section describes both 
techniques. 

Figure 18 shows tt» comfort zone that was derived by Kuznetz wd 
the zone that usually was used in the experiments at SRI. In the 
experiments that are described in the next section the subjects 
universally reported that they were too cold when T^ ^ 16°C or leas and 
the Kuznetz zone was used. The width, slope, and intercept of the 
experimental zone, which can be changed from the keyboard during an 
experiment, were modified so that the subjects reported comfort over the 
c- uplete range of T^ that was encountered in the experimunts. 

A difference in the zones should be expected because of 
different test conditions during Kuznetz* s and SRI 'a experiments. 

Kuznetz used a full space-suit where significant heat was removed by 
airflow and sweat evaporation. At SRI the subjects were enclosed in 
a down insulated garment and there was no air flow; virtually all the 
heat removal was via the LCG. 

The optimum comfort zone is not required to experimentaxly 
verify the method; if the scheme works with one zone, it will work 
with another. Since our Insulated garment did not give the same 
conditions as a space suit we did not feel that an elaborate set of 
tests CO optimize the zone pareuneters was Justified at this time. We 
have the impression that there is no zone that is most comfortable 
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FIGURE 18 COMPARISON OF SRI AND KUZNETZ COMFORT ZONES 


49 




V. '' 


• 9 

J 


C*f ‘‘• 

for everyone, and it aopears simple to experimentally pick the 
parameters to tailor the zone for an individual. This should be a 
subject of further study. 


A. The Iterative Method 

An approximate steady-state AT can be obtained by starting with 

Eq. 13. Assuming that AT* will be near and that the trial T^ will 

equal I_, 
in 

AT* S T^(Cj-l) . (19) 


Substituting in Eq. 14 


41* S 


Where (line 2080) 


„ 0 (n) 

D = m , 

I (n> 

n 

the center of the comfort zone is the line 


( 20 ) 


( 21 ) 


Tc = at* , (22) 

where A^ is the intercept on the vertical axis and Kjj is the slope 
(both are input parauneters; see line 9012).* 

The initial trial T is the T that is a solution to both 
c c 

Eqs. 20 and 22. The algebra to derive the solution is 
straightforward, and the solution is as follows (see lines 9500 
through 9550 and line 9030): 


*In line 9012 and line 9200 K4 = A4 x 9.389. The number 9.389 is 
the slope of Kuznetz's comfort zone and A4 is an input parameter 
to adjust to SRI '3 zone. A5 is a keyboard input parameter for 
adjustment during the experiment (normally set to 0). 
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(23) 


K4 - K4D, 


(24) 




(25) 

I. = 

-B 

■ a--- a 

• - 

q q 

(26) 


2Aq 


The variables I and D are continuously calculated and Eq. 26 gives 
m c 

the initial trial T^. 

The iteration procedure to choose the remaining trial T^ is 

illustrated by the ficticious situation that is depicted in Figure 

19. Based on the present hT and T^ an initial trial T^, T^(1) is 

chosen and the model is used to calculate the first [AT*, AT* (1)1 . 

The point is outside the comfort zone so the next trial [T , T (2)] 

is picked by dropping a vertical line to the zone center and moving 

T (2) down some fraction, H., of the vertical distance, 
c a 

The point *jLT*(2), T (2)1 is not in the comfort zone, so another 

c 

vertical disteuice is measured and T^(3) is moved down Hj. of the 
vertical distance. After the model run, the point [aT*( 3)» T (3)1 is 

found below the zone. The trial T is moved up by half the movement 

c 

fraction that was used to pick T^(3) because there has been movement 

c 

across the zone (each time there is movement across the zone the 

fraction is half the previous fraction) . In our ficticious example 

the fourth trial T^ results in a point in the comfort zone, [aT*(4), 

T (4)1, and this stage of the T selection is complete. 

0 c 

More specifically the rule for the variable movement fraction, 
H^(n), at the nth iteration is (line 9237) 

2Hj.(n-1) if the zone is not crossed 
and 2H^(n-1)< 

if 2H^(n-l) > or s 

.5 H^(n-I) if the zone is crossed (27) 

where H . is some number such that 0 < < 1 . 
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FIGURE 19 A SEQUENCc OF TRIAL POINTS 


52 



The rational for the variable movement fraction is as follows: 
During extensive iterative tests of the scheme it was found that a 
const 2 tnt-movement fraction frequently resulted in Jumps across the 
zone on each new trial with no convergence into the zone. If Jumps 
in the trial T s are sufficiently small, the movement of the point 

G 

will not be enough to cross the zone. However, convergence to the 
zone will be very slow. If the zone is bracketed on each trial and 
the variable fraction is used, the Jump size will be reduced so that 
the point will eventually fall into the zone. When the point remains 
on the same side the fraction is doubled to prevent the fraction from 
being "locked into" a very small value. 

During some series of trials there is not convergence into the 
zone within the allotted number of trials. The reason is unknown, 
but it is suspected that there are situations where the scheme 
described above results in limit>cycles. There was not sufficient 
time to investigate this suspicion, and an investigation of amd 
possible improvement in the scheme should be considered for further 
work. 

B. Single-Pass Calculation 

When there has been no convergence with the multi-pass scheme, 
the single-pass technique is used. It is based on the fact that 
curves of T^ as a function of AT at steady-state have a constant 
slope when the metabolic rate is constant. This can be seen in 
Figure 20. For example, if T^^ is changed from 27®C to 21®C and W = 
500, the change in AT is 0.8°C. Thus the slope in this regior is 
6/0.8 = 7.5. Along the vertical line at W = 500 the curves of 
constant T^^ are equal distance apart, so that the slope is constant. 
The curves when W = 200 are also equal distant but farther apart. 

When W = 200 the slope is greater than when W = 500, 

If T^ is perturbed and the system is allowed to settle (go to 
steady-state) the point (AT, T^^) will move along a straight line if W 
remains constant. During the calculation, two trial T^s, one 
slightly above and one slightly below the present T^, are entered 
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FIGURE 20 SLOPE FOR A CONSTANT METABOLIC RATE. Slope = (27-21 )/0.8 
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Into the model to calculate two AT*s. The slope of the model steady- 

state trajectory is then calculated. Next a line with this slope Is 

projected from the present (AT, Tj^) to intersect with the center of 

the comfort zone. The T^ at the intersection is the new T^. (The 

equation for T^ is in line 9760 where CW is the new T^, SO is the 
c c 

slope, DI is the mean of AT and IM is the mean of T^.) 

Several times in this section the phrase "present (AT, T^)" was 
used. For example, present (AT, T^) is tested against the comfort 
zone. In this test and in the T^ calculations the running averaiges 
of T^ and AT (IM and DM in the listing) are used (not the actual 
present values of T^ and AT) to reduce the effects of noise and 
random variations in these variables. Equation 10 is used to 
calculate the average of T^ and an identical equation with AT as the 
input is used to calculate the average of AT. 

Figure 21 is a schematic of two typical trajectories when new 

T^s are commanded. When (AT, T.) is above the zone the new T will 
c i c 

be less than T, ; when the point is below the zone, the new T will be 

X c 

greater than T^^. When there is a step change in T^ and the 
trajectory is above, the trajectory will swing farther away from the 
zone in the AT direction because AT rapidly increases with the drop 
in T^ and without an initial equal drop in skin temperature. As the 
skin temperature drops, the trajectory will move back toward the zone 
wi^h a decrease in AT. Trajectories below tend to be more parallel 
to the zone. 

Trajectories will remain nearer to the zone (particularly those 
above the zone) if the Tj^ is not abruptly changed to the new T . 

There are a number of obvious schemes for changing T^^ that do not 
involve a step. For example, T^ could be moved along a ramp to the 
new "target" T^. The method that has proven successful in our 
experiments moves T to the target T (the one chosen in the T 
calculations) with the following first order difference equation (see 
the subroutine starting at line 11000): 

T^(n) = (1-A^)T^(n-1)+AJ^ (28) 
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FIGURE 21 TRAJECTORIES OF AT vs Tj AFTER STEP CHANGES IN T^ 
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where is the target and is a constant less than 1 which has 
one of two values depending on whether the trajectory is above or 
below the zone. Equation 28 is Iterated once a second. 

The command temperature will experimentally approach Tj. and the 
time constant of Eq. 28 can be chosen to minimize the distance the 
trajectory of (AT, T^) moves away from the zone. It is our 
subjective opinion that the subjects were more comfortable when this 
scheme was used, and we frequently observed that the trajectory moved 
into the zone and then followed the zone down to the target. 

The next section describes the "closed>loop” experiments where 
the LCG temperature was completely controlled by the computer and the 
subject was exercising on a treadmill. 
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VII LCG AUTOMATIC CONTROLLER TESTS 

A second series of human subject exercise tests was completed 
during the LCG automatic controller progreun. The purposes of these 
tests were to: (1) test the adaptive control program; (2) optimize 
the thermal comfort subroutine; and (3) demonstrate the use of the 
LCG Inlet temperature automatic controller technique during treadmill 
exercise. The computer model actively selected the command LCG Inlet 
temperature during all phases of these tests. 

As was mentioned In the last section, the thermal comfort zone 
as defined by Kuznetz has been found to be Inapplicable for our 
experiments. The subjects tested by Kuznetz were enclosed In a 
sealed suit with ventilation flow that could remove body heat by 
sweat evaporation and convection. Our subjects had no ventilation 
and thus essentially all metabolic heat was removed by a combination 
of respiration and the LCG. Thus, in our case, at any particular 
metabolic rate the LCG must remove more heat for the subjects to be 
comfortable. For a given inlet temperature the AT will be greater 
than given by Kuznetz so that our comfort line appears to hav3 a 
smaller slope. The Kuznetz comfort zone has a Y intercept of '^3*>«6°C 
and a comfort band of ±0.15°C (see Figure 18). The Initial SRI 
comfort zone had a Y intercept of 23°C and a comfort band of ±0.25°C. 
The slope and comfort band of the SRI curve remained constant during 
the various parametric study exercise tests. However, during the 
first series of tests the different subjects tested subjectively 
preferred different Y Intercepts when standing at rest. This 
preference appears to be a function of how well the LCG fits during 
the tests; when the LCG fit a subject tight he preferred that the 
comfort zone Y intercept be Increased as would be expected since the 
thermal conductance between the LCG and the skin Is Increased. 

The thermal comfort subroutine was therefore modified to allow 
the comfort line to be easily varied to accommodate a subject's 
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preference. However, the fact that the comfort line has been found 
to be a function of personal preference and test conditions indicates 
the need for further experimental work. 

Five subject/exorcise tests were completed during the LCG 
automatic controller tests. The results of these experiments are 
presented in Figure 22-26. The same presentation format was used in 
each of these figures. The measured LCG AT (^C) is shown on the X 
axis and the actual LCG inlet temperature (^C) is shown on the Y 
axis. The comfort zone used during each test is presented as the 
solid trace dashed lines in each figure. To provide a uniform 
basis of comparison between these tests and the previously described 
tests the exercise profile given in Figure 14 was used in each case 
studied. Each subject stood at rest for 10 min, walked at 3 mph for 
20 min, and stood at rest for 10 min following the exercise period. 
The thermal conditions during each test are plotted every 100 s in 
Figures 22-26. 

The control program's response was similar during each of the 
five tests conducted. As the subject starts to exercise the LCG 
inlet temperatures versus LCG AT point moves to the right of the 
comfort zone, showing that the subject needs a lower inlet 
temperature to remain comfortable. The adaptive subroutine of the 
control program detects this requirement and altera the control valve 
position accordingly. As the subject continues to work an optimal 
equilibrium point on the comfort zone is reached for each test case. 
The controller automatically adjusts the LCG inlet temperature to 
maintain the desired conditions. Upon cessation of exercise the 
subject's response point falls to the left of the comfort zone. This 
indicates that a higher inlet temperature is needed to maintain 
thermal comfort. Once again this requirement is detected by the 
control program which makes the appropriate valve adjustments to move 
the LCG inlet temperature versus LCG AT point back into the comfort 
zone. 

A second exercise period was Initiated during the test shown in 
Figure 24. As shown In Figure 23 the LCG inlet temperature was 
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LCG DELTA T - 


FIGURE 22 


LCG INLET TEMPERATURE vs LCG AT--SUBJECT SCL/12 






0 1 2 3 4 5 

LCG DELTA T - °C 


FIGURE 23 LCG INLET TEMPERATURE vs LCG AT --SUBJECT PC' /56 
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:rature 



FIGURE 25 LCG INLET TEMPERATURE vs LCG AT--SUBJECT LCL/12 
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FIGURE 26 LCG INLET TEMPERATURE vs LCG JlT--SUBJ£CT GCL/12 
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autoiMtlcally r«a<lju9t«d to • valve position similar to that used 
during the first exercise period. 

A Y Intercept of 23*^C was used during the tests illustrated in 
Figures 22 and 23. The subject tested in Figui'e 23 complained of 
being cold before the start of the exercise. The Y intercept was 
then raised to 2U^C and the test repeated. The repeated test results 
are given in Figure 24. The subject did not feel cold during the 
Figure 24 control period. 

A Y intercept of 26'^C was used during the tests illustrated in 
Figures 25 and 26. The LCG fit these two subjects titter than it 
fit the previous subjects. Initial tests on these two subjects 
showed that they were uncomfortable at rest with lower Y intercept 
values. 



VIII CONaUDING REMARKS AND RECOIfffiNDATIONS 
FOR FURTHER WORK 

A. Concluding Rwaarka 

We have demonstrated that the adaptive-model-reference method is 
a viable technique for the control of LCG temperature. Although the 
algorithms au*e complex, they can be implemented on a reasonably-sized 
microprocessor of a type that probably would be needed for other 
life-support, back pack, house keeping chores. Thus there will be 
liutle increase in the hardware complexity when compared to 
temperature controls that are much less sophisticated. Both the 
house keeping and temperature control should be able to be 
iiq>lemented on a single microprocessor. 

We have developed an experimental microcomputer-driven system 
that can be used to test various algorithms. This system includes a 
through-the-computer loop for the control of the LCG inlet 
temperature. The study developed a servo algorithm that performs 
satisfactorily under the special conditions that are encountered with 
LCG inlet temperature control. We feel that elements of this 
algorithm probably will be included in any operational computerized 
LCG cont4*ol. 

We have developed a simple but accurate adaptive computer model 
that gives outlet temperature as a function of the inlet temperature 
waveform. The model can be run much faster than real-time on a 
microprocessor. It is the heart of the LCG temperature control 
algorithm, and it also might be of value in general computer 
simulations of life support systems. 

There is evidence that the comfort zone varies among 
in viduals. A simple parameter change allows the LCG temperature 
algorithm to be adjusted to virtually any comfort zone. If one life- 
support system were to be worn by more than one subject, a subject 
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Identifying number could be entered to automatically select the 
correct parameters. 

The remainder of this section describes additional work that is 
required or desirable prior to the development of an operational 
adaptive'^odel-reference LCG temperature control. 

B. Recommendations for Further Development 

1. Mechanical Subsystem Improvements 

k nuiid>er of areas were identified during this program where 
additional work is needed to optimize the controller logic before 
proceeding to an actual EVA system. The existing breadboard 
controller can be used for this optimization if several improvements 
are made. 

It was found that the thermal time constant of the system is too 
long. This is primarily caused by the commercial tube-in-shell heat 
exchangers that were used (see Figure 5). These heat exchangers were 
bypassed for the human tests discussed previously. However, this was 
only nmu*ginally effective due to the large mass of water in the 
system. These heat exchangers should be replaced by aerospace 
quality compact llquld-to-liquld heat exchangers that are available 
from a number of sources. 

The proportioning valve must also be replaced. The present 
valve, which was borrowea from NASA, leaks both internally and 
externally, responds slowly to conneunds to change position, and has 
excess mechanical hysteresis. 

If the above hardware improvements to the breadboard system are 
implemented, tue following test series, with accompanying analysis 
auid control logic optimization, should be performed. 

2. Tests Using an Insulated Overgarment 

Both the steady-state and transient controller responses should 
be investigated further. Limited human evaluations were performed 
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during the current program. In particular, very few repetitive runs 
of the type needed tc develop statistically significant data were 
performed. It is expected that the controller performance can be 
improved significantly by these tests. 

3. Whole Body Heat and Mass Balance 

Following the tests discussed above, a series of human 

experiments should be performed to completely quantify all aspects of 

the controller's physiologic performance. These tests would be 

q 

similar to those performed by Webbon wherein subjects were enclosed 
in a sealed environment simulating a space suit. The subject's 
thermal state would be continuously assessed by measuring the suit 
air and water inlet and outl.^c conditions, as well as physiological 
parameters such as heart rate, actual metabolic rate, and various 
body temperatures. The subject's precise thermal state can be 
calculated from these data. Repetitive runs would be performed using 
a variable metabolic rate to completely characterize the controller's 
performance. 

4. Model Coefficient Refinement Tests 

A preliminary model coefficient parametric study was performed 
as part of this investigation to determine the number of taps, the 
adaptive coefficient, and the mean coefficient values that should be 
used in the control program. The procedure and the results obtained 
during this parametric study are described in Section V.B. of this 
report. As discussed earlier, the parameter values used led to an 
average error of about 0.4°C between the predicted and measured LCG 
AT values for a series of four exercise tests. It is recommended 
that additional parametric calculations be made to verify and refine 
the previous coefficients. 

The measured LCG AT values obtained during the insulated 
overgarment tests will then be used in a parametric study to 
determine the coefficient values that provide the best calculated LCG 
AT values. The following parameter values will be studied during the 
parametric portion of this test series; 
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• Number of taps: 2 , 4, 6, 8 

• Adaptive coef f io: int : 0.15, 0.10, 0.075, 0.05, 0.025 

• Mean coefficient: 0.1, 0.2, 0.3. 

The optimal set of these coefficient values will then be incorporated 
into the automatic controller prograun. 

5. Dither Signal Optimization Tests 

An additional exercise test series should be performed to 
optimize the amplitude and frequency of the dither signal required 
for the adaptive control model. A random dither signal of ±1°C, with 
a dither interval greater than 30 s was used in all previous tests. 
System performance and subject thermal sensitivity will be determined 
with the following dither signals: 

Amplitude(°C) Interval (s) 


0.25 

20 

0.25 

60 

0.25 

120 

0.50 

20 

0.50 

60 

0.50 

120 

1.00 

20 

1.00 

60 

1.00 

120 


Two subjects will perform the metabolic profile shown in 
Figure 27 during treadmill exercise with the above dither signals. 

The LCG inlet temperature will be under full automatic control during 
these tests. Metabolic rate will not be measured directly but will 
be determined using standard correlations with the known treadmill 
speed. The LCG will be an Apollo-type garment. The subjects will 
wear an insulated down suit over the LCG. The subject’s thermal 
state will be determined using the same thermal -comfort index used by 
Kuznetz. No other physiological data will be measured at this time. 
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Control system performance will be assessed using the LCG inlet 
and outlet temperature profiles during the various tests. A dither 
amplitude and frequency will be selected which optimizes the system 
control characteristics and minimizes the subject's sensation of 
thermal change. 

6. Hot-Cold hater Inlet Temperature Control 

Use of optimum control techniques and a more rapid return to the 
comfort zone would be possible if a source of hot water (in addition 
to the cold water from the heat exchanger) were available. Then 
could be rapidly driven in either direction and the overshoot could 
be minimized. Optimum control could minimize the time and distance 
away from the control zone. 

There should be a study of hot water sources using heat from the 
LCG. For exanq>le, water from the LCG irtien it is operating at the 
upper range of T^ could be stored in a hot-water reservoir, and water 
from the heat exchanger could be stored in a cold-water reservoir. 
Algorithms for deciding when to withdraw, the rates for withdrawal 
from either reservoir, and the replenishment rates with typical 
metabolic rate profiles should be developed. Reservoir sizes should 
be studied along with control algorithms that utilize a hot/cold 
water control. 

7. Algorithm Improvements 

The control algorithms that are described in this report are 
preliminary developments and additional refinement will be required 
prior to any operational use of these techniques. For example, the 
convergence in the T^ calculations should be better understood amd 
perhaps a modified technique should be developed. It is important to 
conduct an extensive set of experiments (similar to the proposed 
model experiments) ti optimize the parameters in the T^ calculation 
and in the computer servocontrol loop. 
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There au?e over 20 input variables for the servo and 
calculation algorithms. Our data base was too small, and we did not 
have time for more than cursory optimization. 

8. LOG Transient Responses 

Kuznetz's work and the T^ control algorithms that were developed 
during the present study are based on a constant T^, an LOG at 
steady-state, and a constant metabolic rate. None of these are 
strictly true in practice. Kuznetz's curves and comfort zone assume 
steady-state, and the comfort region when the metabolic rates are 
changing is not known. What are the most comfortable T^^ versus AT 
trajectories? Eventually Kuznetz's work must be extended to include 
a dynamically changing environment. 
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Appendix 


PROGRAM LISTING 



The following is a listing of the Microsoft BASIC version of the 
program that is used to control the LOG temperature. This program is 
compiled (Microsoft BASIC coiqpiler) for the experimental studies. 

The parameter values can be obtained by comparing lines 30 and 

34,000. 
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page t3 


5 CLEAR :;OOJDIh CT • 1 0 0 > , In ( 1 00 ) , Dh ■! 1 0 0 > 

7 AF=n ; AD=i !L3:i-=o ? A3:i=o ; LT»=" 0 " f jo::=o : j^::*o !L7:w=o JFM 


0 } F5%=o ! J8;:=n ;f9:.*o : bx%=o 


10 CL5 

IS OH ERROR GOTO 900 

20 PRIWT-THIS PROGRAM DOES THE CLJSED-LOOP COHTROL f PRINT 
IS 1HPUT"PPESS XNTEP' TO CONTINUE" JA.'CLS 

30 READ CU , BU » BO , C1> > C 1 « , DS , TD': , PH > R » AL , AN > LOi: , L3>: » BN , BL , AS , HC > Si;;! , S J% , HLa , HO « AB» 
»CC,A6.BC1',A7,A8,BA»>L8:v>A0,A<4,AA»BG$>BI«BJ,BH$,AX,AY,AH 

3S IC=AA/ ( 9 . 38<»*A'» ) J AM*AN»L3Z } BM»BN*L3X 1 A1 *A4»9 . 389* < IC-BL ) t A2*A^»9 . 389» ( IC+BL ) : 
A®*A0+m5 

■^0 INPUT"ARE all the default UALUES to be USED (Y/r!)"5Af ;CLS;PRINT;1F A»»"N" the 
N ag*o;gosub i30oo;else ag»i 

“IS FOP N %=1 TO ^SREAD CP(N-;);NEXT NrilFOR N =1 TO R J C ( N ) *0 5 NEXT N?IF R=4 AND AB»*" 

Y" THEN FOR N=1 TO ^IRFAD C(N' :CZ':N>»C(N) JNEXT N 

SO GosuB iiooo;gosub 8000:GOSUB 6000 :GOSUB 7ooo:gosub sooo 

S5 jo:i*Jo:>i:iF jo%; *lo% then jo%*o 

57 J3:i=J37:+l :IF J3:1=L3'. then J3%=0 

eO IF AF=0 THEN PPIK'T 8 321 , "CT" J : PRINT TAE<12>"IT"J TAB<2^ ) "OT" ? TAB ( 36 ) "HT" }TAB 
/i»8)"F" 

70 IF AF»0 THEN PRINT 8 385, CTJ .‘PRINT TAEUDIT? TAB(2^)0T;TAB(36)HT;TAB(48)F 

80 IF FL%=1 AND AF*0 THEN PRINT" Jl*" } JIX 

100 IF A3%=0 THEN COSUB 1000 

no IF LT#»"0" THEN 500 

130 IF JOXsl THEN GOSUB 3000 (GOSUB 2000 

1^0 IF J0Xx2 AND AFsl THEN GOSUB 14000 

145 IF J0%a3 AND FLX=1 THEN GOSUB 12000 

150 IF <J0X=4> OR (FK*>1 AND J0X<L0Z-1) THEN GOSUB 9000 

500 A»*"“ 


510 A1‘*INKEY» 

520 IF Af»"C" THEN AD'.= 1 J GOSUB 13000 
530 r«* TfsTIMEt THEN 500 
540 T.=TIME»:COTO 50 
550 POKE 16425,1 

560 CLS:iNP'JT"RECORD THE INLET TEMPERATURES ^Y/N)"JA<SIF A»»"N" THEN 640 
570 PPINT:INPUT"ENTER the file NAME";NF4 
580 OPEN"0",l,NF»+":i" 

590 PRINT*1,J1X 
600 FOR NX»1 TO JIX 

610 prtnt*i,ct(nx) ;im<n:o jdm(nx) 

620 i T NX 

630 close:jix*o 

640 PRINT :iNPUT"CONTINUE THE RUN ( Y/N) " } AU* J IF AU»»"Y" THEN 50 
650 END 

900 RESUME 910 

910 JlX=JlX-i:iF OT-HT=0 THEN F=F+BS*(CT-IT) :FS»F;IF IT-CT<0 THEN FH*1 

920 IF IT-CT>0 AND OT-HT-0 THEN FH=2 

930 IF OT-HT»0 THEN 7050 ELSE GOTO 550 

10 00 REM 5 THIS SUB INIaLIZES IM AND OM 

1010 L7X-=L7%+i;iF L7X = 10 THEN A3X = 1 

1(20 IM*IM+IT:OM»OM-t-OT 

1030 IF A3X=0 THEN 1050 

1040 IM*IM/10:OM*OM/10 

1050 RETURN 
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ORIGINAL PACE 13 
OF POOR QUALITY 


I»100 RErttTHIS SUB CALCULATES THE COEFFICENTS 

Z0Z7 IF (DT~DE>C2: 1 THEM KC=SGp! ( DT-DE > ELSE KC=DT-DE 

2030 FOP NC=i TO P 

2050 C(HC)=C(NC>+AL»KC»XC<NC) 

2070 NEXT HC 
2080 DC=OM/IM 
2090 RETURN 

3000 REM : THIS SUB DOES THE FILTERING 
3010 FOR Nr==R TO 2 STEP -1 
3025 XC^NF)=XC(N^ 1) 

3030 NEXT NF 
3035 XC<1)=H1 
30^0 OE=0 

3050 FOR NF= 1 TO R:DE=DE+XC <NF)»C (NF> : NEXT NF 
3C60 DE=DE-^<DC*1)kIM 
30*0 RETURN 

5000 REM: THIS SUB-PROGRAM SETS THE VALUE POSITION 
5010 REM: variable 'PV' is iHl VALVE POSI'^^ION 

5012 iP=iP+i:iF ip :ph the^» 5120 else if ct<: rt then ip=o:rt=ct 
5020 PA=223 

5030 MS=INT<PV/256 + *00 
50^0 OUT PA*16*l-»-nS 
5050 OUT PA,f1S 

5060 MI=INT( ^PV/16>-163KiiS+.001) 

5070 OUT PAtl6»2+MI 
5080 OL^ PA, MI 

5090 LS=INT(PV-256*M l6*MI+*00l) 

5100 OuT PA,16»3+LS 
5110 OUT PA,LS 
5120 RETURN 

6000 RE*.: THIS SUB READS THE TEMPERATURES 

6010 OUT 223,128 

6015 FOR MT=3 TO 10:NFXT MT 

6020 FOR NT^=0 TO 4 

6022 IF FG%=1 AND NT%=^ THEN 6050 

6025 FOP MT=1 TO 10:NEXT MT 

6030 OUT 223, NT/: 

6032 FOf MT=l lO lOZNEXT MT 
6035 OUT 223, 12P 
6037 FOR MT=1 TO lOINEXT MT 
60^0 TT<NTr:^*l)=INP(223) 

6045 IF NT% :4 AND TT ' NT%-M ) ;= 1 15 THEN TT(NT%+1 )=VT(NTX^-1 > 

6047 vt(nt::+i)=tt<nt%+i) 

6050 NEXT NTX 

60J5 IF TT(5>: lOO “tHEN FG%=1 
6060 HT='TT(3)-51*4«^)/1*7P 
6080 IT=(TT(2>-4P*95)/1*82 

6000 0T=(*T(1)-51*49)/1*79 

6001 IF FG>:=1 THEN DT=OT-‘ T ELSE DT=(TT(5>-52. ( ( (0T+IT)/2)-8) )/22*48 

6093 IS=ITtOS=OT:iT=(ir+OT-DT>/2:iF IT-: 0 THEN IT*IS 

6006 0T=^IS-»‘0T+0T)/2:Tr OT !2 THEN 0"'=OS 

6100 IF fg::=i and dt . ^ then fc::=o 

6105 1FJ3%^Q THEN Dh^ . " 1-A. ))tDi : IM=AM«IT+< 1-AM ) wIM :0M=AM»^*T+ ( 1 -AM) »0M :H1=IT 

-TE-f l-BM;*Hi:iE=IT 
61 RLTURN 
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ORIGINAL PAGE 13 
OF POOR QUALITY 


■’00 0 REM : THIS SUE GEHEPATES THE MALVE POSITION. 

~010 IF CT=ST THEN 7033 ELSE ST=CT 
'015 IF C1*="N" AND CT-IT>=0 THEN 7033 

7020 F=F+(CT-IT)/(OT-HT) :FS=F:IF IT-CT.O THEN FH=1 ELSE IF IT-CT>0 THEN FH 
■^030 GOTO 7950 

7033 IF CT-IT>=0 THEN BS=BU:G0T0 7 0'90 
7036 GOTO 7200 
7040 F=F^eS»(C IT) 

7050 PM=0 

7053 IF FH=1 AND IT-CT? =0 THEN F=FS:FH=0 

7056 IF FH=2 AND IT-CT<=0 THEN F=FStFH=0 

~060 IF FI THEN F=1 

•^070 IF F' 0 THEN F=0 

~080 FOR NS5:=1 TO 3 

'090 PV=PM+(FC(4-NS2) )*CP<NS%) 

7100 NEXT NSZ 

7110 pv=pv*cp<4);pv=int(p*.»);if pv<i then PV=1 
7115 IF PV>4095 THEN P')=4095 
'120 RETURN 
7200 BS=BD 

7210 IF C*="N** THEN 7040 

7220 IF (IT-CTXl THEN 7040 

7230 F*F-BS:G0T0 7050 

8000 PEMITHI SUB PICK THE DITHEP 

8010 IF DS=0 THEN CT=CU;GOTO 8070 

8020 JD%=JD!2+i:iF JDX=TDX THEN JDZ=0 

8030 Tl4=MID4(TT,14»l):iF T14="7” THE< CT=CU+DS/2 IGOTO 8070 

R04C IF Tlt="8" THEN CT=CU:G0T0 8070 

6050 IF T14=”9" THEN CT= CU-DS/2:G0T0 8070 

6060 IF JD%=1 THEN XD=RNO<2) :CT=CU+DS» ( <XD-1 )-.5) :CM=CU 

8065 IF CWOCU THEN CT=CU+DS»< (XD-l )-.5) :CN=CU 

8070 RETURN 

8500 pen: THIS SUB CHANCES THE DITHER SIZE 
8510 CLS: INPUT-ENTER THE DITHER SIZE" ;DS1 PRINT 

8515 INPUT "ENTER THE TIME IN SECONDS BETWEEN DITHER CHANGES" ;TDZ:CLS 
8520 RETURN 
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0R13!NAL PAGE 
Of POOR QUALITY 


«>00 0 REM : THIS SUE PICKS THE CT DURING CLOSED-LOOP OPERATION 
^OOS IF LTt="A“ AND FM%=0 THFW 9260 
®010 IF FM%=1 THEN ®060 

<>012 IF im: aa-a-^»o*38o*dm+a5 THEN sh%=sj:l ELSE sh::=si% 

9015 J4%=J^:i:-H:iF J^Z > 'riZ and F9Z=0 then 9260 

9020 IF IM>A1-A^*9.389r .‘M+AS AND IM<A^-A^*9 . 389*DM+A5 THEN 9260 

9025 IF F9Z>0 AND ^SHZ THEN 9950 

9027 IF BG$=”y THEN 960i 

9030 D0=DC:G0SU8 9500:PPINT 8 63,”»" 

90^0 FOR NHZ=i TO RtyH(NHZ>=xc<NHZ) :CH(NHZ)=C(NHZ > :next nhz:dd=dc:ih=im:h2=hi:h 2 
=cu-iT+ • 1 -BM ) mh 2 : di=dm : hf=ho : cy=cu 

9050 FOR NHZ=1 to HLZ 

9060 IF TSOTIMEf THEN FMZ=1:G0T0 9260 
9C70 FOP NIZ=P T02 STEP -1 
9080 XH(NIZ>=XH(NIZ-1) 

9085 NEXT NIZ 
90®0 XH(1)*H2 


9100 DF=0 

9110 FOR NIZ=1 TO R:0F-DF+XH(NIZ)*CH(NHZ) tNEXT NIZ 
9112 H2=(1-BM)«H2 


911^ NEXT NHZ 

9115 CD*DD»(l+HC»(IH-CU)):iF CD<1 THEN CDsCC 

9120 DF*DF+(CD-l)»CU:iF (DI-DF)C2>2*25 OR DF<0 OR DF>5.5 THEN DF=(CO-))*CV 
9150 IF F5Z*1 THEN 9250 

9200 IF LCZ=0 AND CM>=AA-A^*?.389*DF^A5-A6 AND CV <A2-A4*9.389*DF+A5 THEN IF 
’•Y*' GOTO 9600 ELSE GOTO 9250 

9210 IF LCZ=1 AND CU<=AA-A^»9*389)rDF^A5+A6 AND CV>Al-A4»9*389xDF+A5 THEN IF 
••r* GOTO 9600 ELSE GOTO 9250 

9215 IF BC$=“Y** AND CV>AA-A^*9,839*(DF+A7>+A5 AND CU< AA-A^»9.389»OF-A7)+A5 
IH-CV A8 THEN 9600 

9220 CV=CU+HF)r ( AA-A4»9 . 389*DF^^A5-CV ) 

9225 LDZ=LCZ 


BC$= 


BC$= 


AND 


9230 IF CV<AA-%^»9»389»DF-^A5 THEN LCZ=0 ELSE LCZ=1 
9235 IF LDZ :>LCZ THEN HF=HF».5 

9237 IF LDZ=LCZ THEN HF=HF*2:iF HF>HD THEN HF=HD 

9238 J8Z=J8Z+i: IF J»Z>*L8Z THEN 9600 


02^0 GOTO 9050 


9250 cr»cu:bxz=o:f9z*o:j8z=o:f5z=o:fmz=o:cu=cu:print b 63^** *':x3=df/.05+i:y 3=^6- 
^CU-10>»2*^:iF X3<0 THEN X3=0 ELSE IF X3>127 THEN X3=127 

9251 IF Y3 0 THEN Y3=0 ELSE IF Y3>-<»6 Y3=^6 

9252 IF AF=1 AND FR=0 THEN RESET (X^rY4) 

925^ IF AF=1 AND FQ=0 THEN RESET(X^+1 , Y^) 

9256 IF P0INT(X3,Y3) THEN FP*1 EL" E FR=0 
9256 IF POINT (X3+lrY3) THEN FO=l Ei.SE FQ=0 

9259 X^=X3;Y^=Y3:IF AF=1 then SET(X3,Y3) :SET(X3^-1 »Y3) 

9260 RETURN 
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cr.il:'NAL ^3 

Of POOR QUALITY 


0500 A0=A^»9.389*00»HC 

<>510 BQ=A^»o , 389« ( l-DD-D03rHC*IH > -1 

9520 CQ=AA^A5 

9530 D0=BQ»BQ-4*AQ»CQ 

95^0 IF DQ<0 THEN 9560 

9550 CV*(-BQ-SQR(DQ>)/(2*AQ):G0T0 9570 

9560 CV=HO*(AA-A^»9.389»DM'^A5> 

9570 IF CV^^5 THEN CV=^5 

9580 J^Z=0:iF CV-IH>0 THEN LCZ«0 ELSE LCZ=1 
0590 RETURN 

9600 IF BA»=*'N» THEN CM=CU:F5Z=i: GOTO 90^0 
0605 PRINT 8 63. 'V: FOR NJZ=1 TO 3 STEP 2 
9610 FOP NHZ=1 TO RtXHf NHZ)=XC(NHZ) INEXT NHZ 
<^630 H2=H 1 : H2=N JZ-2+ ^ 1 -BM ) »H2 : OI=DM t HF=HD 
96^0 FOR NHZ=1 TO HL2 
9650 FOR NIZ=R TO 2 STEP-1 
9660 yH(NIZ)=XH(NIZ-l) 

0670 NEXT NIZ 
9680 XH(1)=H2 
9690 OG=0 

9700 FOR NIZ=1 TO R :OG=DG+XH( NIZ )*CH( NHZ ) t NEXT NIZ 
9710 H2*a-BM)«H2 
9720 NEXT NHZ 

0730 CG=DC*<l-^HC»(2-NJZ)):iF CG<1 THEN CG*CC 
07^0 OG ( N JZ ) *DG+ ( CG-1 ) » ( IH+N JZ-2 ) 

9750 SD=2/(DG'3)-DGa)) 

975^ NEXT NJZ 

0760 CW= ( SD» AA-A^»o , 309* ( SD*DI-IH ) +A5*SD ) / ( SD+A^*9 ♦ 389 ) 

9765 IF CV>CU KH=KL ELSE KH=KM 

0770 IF (SD<-A^«9,389 and J8Z>=L8Z) OR (BGO=”Y“ AND SD<-A^*9.380) THEN CV=CW:DF 
(AA+A5-CM)/(A^»9.389> :ELSE IF SD! =-A^*o,389 AND J8Z>=L8Z THEN F9Z=i:G0T0 9900 
0780 IF CV:Ao then CV=A9:F5Z=1 :G0T0 90^0 
0785 IF CV: AH THEN CV=AH:F5Z=1 tGOTO 9040 
9790 GOTO 0250 

9900 IF BA$="N*‘ THEN CU=CU:F5Z=1 tGOTO 9040 

9910 SE=-BI*A4»9 . 389 1 CV= ( SE*AA-A4»o , 389* ( SE*DI-IH ) +A5*SE ) / ( SE-^A4*9 . 389 ) : F9%=2 
9920 F5Z=1:G0T0 9040 

0050 IF CU<AA-A4*Q*380*Dh'«'A5 THEN CX=-1 ELSE CX=+1 

9955 CU=CU+BJ*SE*CX:F0%=1 

9960 IF CU <A9 THEN CU=A9 

0965 IF CU>AH THEN CU=AH 

0970 GOTO 9260 

11000 REM : THIS SUB DOES THE CT FILTERING 

11010 IF BH^=“N** OR F9Z=1 OR LTt- >”C» THEN 11040 

11020 IF BXZ=0 THEN CU=CY:BXZ=1 

11025 IF CU CR then AT=AX ELSE AT=AY 

11030 CU^d-ADxCU-^'ATxCR 

11040 RETURN 

121-00 REM : THIS SUB STORES VARIABLES FOR RECORDING 
12010 ji:;=jiz-d 

12015 IF AF=1 THEN PRINT 81 0 1 , " J1 = ** ? JIZ 
12020 CT( Ji::)=CT;iMMi;:>=iMrDM( JiZ)=DM 
12030 RETURN 
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Cr ' 

Of 


13000 

REM : THIS : 

BUB Cl 

13010 

CLS 

:: PRINT 

8 o2 

130^0 

At = 

. M «: 


13050 

Af = 

INKEYT 


13052 

IF 

A^=”b” 

THEN 

IS SHIFT 

B 


13055 

IF 

A^=”0” 

THEN 

0 13500 



13060 

IF 

A1^=”T” 

THEN 

13062 

IF 

A^=”r” 

THEN 

:bi 




13H65 

IF 

A^-”8” 

THEN 


THEN INPUT”IS THE SLOPE PREDICTION MODE TO BE USED” ;BG$ :RE«- <Br 


THEN INPUT-ENTER THE DT INTERCEPT VALUE” : IC :A4»AA/ 38<>»IC >: GOT 


THEN PPINT-THE VALUE OF A4 NON IS” J A4 : INPUT”ENTEP THE NEU SLOPE 
CHANGE-FACTOR. A^ CENTER 0 IF THERE IS TO 

BE NO CHANGE' ”:AZ: IF A2 0 THEN A^=AZ:G0SUB 13500:ELSE GOSUP 13500 
IF At=”G* 


13070 
130*"2 IF A^-”j* 
IF AT*”*^* 
IF A^=”D* 


THEN FLX=1 

THEN INPUT-ENTER THE SIZE OF THE CT-SERVO STEP”;BJ 
THEN INPUT-ENTER THE VERTICAL INTERCEPT ^ AA”;AA:G0SUB 13500 
THEN GOSUB 8500 

IF A^=”h* THEN INPUT-IS THE CT-FILTER MODE TO BE USED (Y/N)”;BH$ 

13085 IF At=”6” THEN INPUT"ENTER THE LOMEST CT”;A9 

13090 IF A^*-C- THEN INPUT”ENTER THE NEW UP SERVO GAIN CONSTANT” ;BU: PRINT: INPUT* 

ENTER THE NEW DOWN SERVO GAIN CONSTANT” JBD 

130«>5 IF Af=”5” THEN INPUT”IS CT TO BE ADAPTED WITH OUT THE MODEL BEING ADAPTED 


130*^5 

13080 

13082 


(Y/N)- 

!BAf 


13100 

IF 

At=”E" 

THEN 

13102 

IF 

AS=-a" 

THEN 

13105 

IF 

AS--4” 

THEN 

(Y/N) 



13110 

IF 

A^=”F” 

THEN 

)-:cif 



13112 

IF 

A$="l<” 

THEN 

13115 

IF 

At^=-3” 

THEN 

EN 




IH-CV^ 

A8) 

> 

CD 


13120 

IF 

A^*-H- 

THEN 

13122 

IF 

AS*”i” 

THEN 


ABOVE THE 
COMFORT-ZONE” ;SJZ 
13125 IF A^*”2” THEN 
(A7)-:A7 

13130 IF A^=”X” THEN 
13132 IF At*”q” THEN 
13135 IF At=”l” THEN 
131^0 IF Af*”A” THEN 
131^2 IF At = ”fi” THEN 

I :low the 
COMFORT-ZONE" ;si!i: 
131^5 IF At=”0” THEN 
13150 IF A$="M” THEN 
13152 IF A^=”t” THEN 
13155 IF A^=:”0” THEN 

”:ap$:if ab$^”Y" and 
n r:c<nz)=o:next nz 


550 

INPUT-ENTER THE DOWN CT-FILTERINC CONSTANT”: AX 

INPUT-ARE THE EXTRA LOOK-AHEAD STOPPING CRITERIA TO BE USED 

INPUT-IS THE INITIAL F TO BE USED AFTER A DOWN COMMAND (Y/N 

INPUT-ENTER THE UP CT-FILTERING CONSTANT”; AY 

INPUT**ENTER THE VALUE ^F A8> <THE SECOND LIMIT OPERATES WH 


INPUT**ENTER THE NEW SERVO HOLD TIME”: PH 

INPUT-ENTER THE NUMBER OF STEPS TO HOLD CT WHEN (DM,IM) IS 


INPUT-ENTER THE VALUE OF THE CONSTANT FOP THE SECOND LIMIT 

INPUT-IS THE STEPPED-COWN COMMAND TO BE USED ( Y/N )” ;C$ 

PRINT” J^X=” : : INPUT-ENTER : J^X 

INPUT-ENTER THE NEW VALUE OF THE ADAPTIVE CURVE OFF-SET” ;A5 
INPUT-ENTER THE VALUE OF THE ADAPTIVE CONSTANT” ;AL 
INPUT-ENTER THE NUMBER OF STEPS TO HOLD CT WHEN (DM,IM) IS 


INPUT-ENTER THE NEW DEAD-ZONE C0NSTANT”:A6 
INPUT-ENTER THE VALUE OF THE MEAN CONSTANT": AM 
INPUT-ENTER THE UPPER CT-LIMIT-:AW 

INPUT- ARE NON-ZERO INITIAL VALUES FOR C(N) TO BE USED (Y/N> 
R=^ THEN FOR NZ=1 TO ^ : C ( NZ > =CZ ^ NZ ) : NEXT NZ:ELSE FOR NZ=1 T 
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»- V - 


13160 IF A1-="P" THEN INPUT"ENTER THE VALUE OF THE DC PROJECTION CONSTANT” tHC 

13165 IF AO=”N" THEN INPUT”ENTER THE NUMBER OF CALCLATIONS FOR LOOK-AHEAD” JHLX 

13170 IF A»=”W" THEN INPUT”ENTER TYPE OF DISPLAY (1 OR 2) ” J AFtAF*AF-l 

13175 IF A»*”I” THEN INPUT“ENTER THE NUMBER OF STEPS TO HOLD A NEW CT WHEN (DM, 

IM) IS 

BELOW THE COMFORT-ZONE” J Sit 

13180 IF A»*”L” THEN INPUT”ENTER <C> FOR CLOSED-LOOP OR <0> FOR OPEN-LOOP CONiRO 
L OR <A> 

FOP OPEN-LOOP CONTROL WITH FILTER ADAPTION” JLT» 

13185 IF A*="V” THEN INPUT”ENTER THE CT ADAPTION CONSTANT” 5HD 

13190 IF A»=”B” THEN INPUT"ENTER THE NUMBER OF SECONDS BETWEEN FILTER CALCULATIO 

NS"?L0X 

13195 IF A$*”U” THEN INPUT”ENTER THE HI-PASS CONSTANT IN 1/SEC. ”; BN :BM=BN»L3Z 

13200 IF A»s”R” THEN INPUT”ENTER THE NUMBER OF TAPS”5R 

13205 IF A**”Z” THEN INPUT“ENTER THE HI-PASS CONSTANT” ;BM 

13210 IF A*=”K”THEN INPUT”ENTER THE DELTA T LIMIT” ;BL:A1»A4*9.389*(IC-BL) lAZ-A-^K 
9.389*(IC+BL) :IF BL<A7 THEN A7*BL/2 

13215 IF A$»"Y” THEN INPUT”ENTER THE MEAN CONSTANT IN 1/SEC.”;AN:AMsANxL3X 
13220 IF A»*”S” THEN FL2=0 

13225 IF A**”J” THEN INPUT”ENTER THE SECONDS BETWEEN MEAN CALCULATIONS” JLSZIINPU 

T“CHANGE AM AND BM (Y/N) ” JAACIF AAi*”Y” THEN AM*ANxL3%:BM»BNxL3Z 

13230 IF AG*0 AND A»<>"0” THEN 13040 

13240 IF A»=”” THEN 13040 

13250 CLS : AG* i: RETURN 

13500 IC=AA/<A4*9.389) 1 A1*A4*9,389*(IC-BL) :A2*A4*9.389*(IC+BL> 

13510 RETURN 
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1^000 REM : THIS SUB PLOTS THE CRT 
1^010 IF POINT< 120,^6) THEN 1^250 
1^020 CLS 

1^030 FOR MCZ=10 TO 3^ STEP 12 
1^040 SET(0rMC%) 

NEXT Mc:: 

14050 SET (0,47) 

14060 FOR MCZ=1 TO 120 STEP 10 
14070 SET(MCZ,47) 

14080 NEXT MC2 
14090 FOR MC:i:=0 TO 46 
14100 SET(1,MC%) 

14110 NEXT MC% 

14120 FOP MC%=0 TO 120 
14130 SET(MC%,46) 

14140 NEXT MCZ 
14150 FOR MC2:=1 TO 120 
14160 MD=(MC%-1)**05 

14170 Y 1 =46- ( A 1 - A4»9 • 389*MD+ A5- 1 0 ) *2 ♦ 4 
14180 I' YKO OR Yl>46 THEN 14200 
14190 S (MCZ,Y1) 

14200 Y ^=46- ( A2-A4*9 * 389*MD+A5-1 0 ) *2 * 4 
14210 IF Y2<0 OR Y2>46 THEN 14230 
14220 SET(MCZ,Y2) 

14230 NEXT MCZ 

14240 FET(X4,Y4>:SET(X4+1,Y4) 

14250 ^ESZT(X5,Y5) 

14260 X5=DM/*05+1:Y5=46-(IM-10)3ic2*4:IF X5<0 THEN X5=0 ELSE IF X5M20 THEN X5=120 
14265 IF Y5<0 THEN Y5*0 ELSE IF Y5>46 THEN Y5=46 
14270 SET(X5,Y5) 

14280 PRINT 037,INT(CU*lOO)/lOOr* “ ;INT<IM*100 )/100 J ” ** }INT(DM*100 )/100 

14290 RETURN 

34000 DATA 21 , ♦ 075, . 025, ”Y” , ''Y” ,2 ,60 ,5,4 , . 05, ♦ 01 20 , 1 , * 005, *25, 0 , * 017, 15, 15,6, *5 
,“Y’M,08, .5,”Y”, .15,4, ••Y’M 0,12, .35,23, "N'M. 5, .6, "N’* , .004, .02,21 
35000 DATA -199.25,-802.063,4912.55,34.6582 
35010 DATA -.491 ,. 0089, . 115, . 187 
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